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ABSTRACT 


The  Army  RAH-66  Comanche  Helicopter  is  currently  undergoing  developmental  flight 
testing.  The  empennage  of  the  aircraft  is  experiencing  buffeting  where  the  horizontal  and 
vertical  tail  vibrate  at  resonant  frequencies.  These  high  buffet  loads  are  manifested  in 
higher  than  anticipated  fitting  loads,  particularly  on  the  tail,  and  vibrations  in  the  crew 
stations  and  at  the  nose  cone  where  the  targeting  sensors  are  located.  Significant  effort 
has  been  devoted  to  identifying  the  sources  of  excitation  and  the  nature  of  the  structural 
response.  This  thesis  determines  the  location  and  magnitude  of  empennage  vibratory 
airloads.  Because  the  nature  of  the  excitation  is  a  random  function,  spectral  analysis  is 
used.  To  obtain  the  loads,  a  three  step  process  was  utilized.  First,  from  aircraft 
differential  pressure  transducers  and  accelerometers,  the  spectral  content  of  the  response 
and  excitation  was  determined  Then,  using  a  NASTRAN  model  modified  to  replicate  the 
flight  test  aircraft,  frequency  response  functions  were  determined  between  selected  points 
on  the  aircraft's  tail  and  the  accelerometers.  Finally,  using  this  information,  a  solution  was 
obtained  for  the  vibratory  airloads.  Having  provided  information  on  the  nature  of  the 
driving  forces,  structural  modifications  can  be  made  that  move  the  natural  frequencies 
away  from  the  frequencies  of  the  applied  airloads. 
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I. 


INTRODUCTION 


A.  GENERAL 

The  RAH-66  Comanche  is  the  United  States  Army's  newest  light/attack  armed 
reconnaissance  helicopter.  It  is  designed  to  be  a  combat  multiplier  on  the  joint  and 
combined  arms  battlefields  of  the  21st  century.  The  aircraft  will  be  fielded  to  both  corps 
and  division  cavalry  and  attack  battalions/squadrons  with  the  assigned  missions  of  armed 
reconnaissance  and  attack,  both  with  an  embedded  air  to  air  capability.  Mission  profiles 
include  day,  night,  reduced  visibility,  from  ships,  and  in  desert  and  non-desert 
environments.  The  aircraft  is  shown  in  Figure  1.1. 

Comanche  will  contribute  to  combat  effectiveness  by  its  ability  to  fight  for 
information  and  disrupt  the  enemy's  tactical  intelligence  collection  effort  in  the  counter 
reconnaissance  role.  Key  attributes  of  the  aircraft  include  a  five-blade  bearingless  main 
rotor ,  F  ANT  AIL  anti-torque  system,  composite  materials,  Low  Observables  (LO),  an 
advanced  pilotage  and  targeting  system,  and  an  integrated  digital  mission  equipment 
package  (MEP).  These  digital  systems  will  allow  it  to  serve  as  a  forward  data  fusion 
center  and  provide  near  real  time  information  to  commanders  at  all  levels. 

Designed  for  two  crew  members,  it  is  capable  of  being  operated  by  a  single  crew 
member  from  either  station.  During  missions,  the  pilot-vehicle-interface  (PVI)  is  designed 
to  support  low  workload,  high  situational  awareness,  human  error  resistance,  and 
maximum  mission  productivity.  [Ref.  1] 
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Figure  1.1:  US  Army  RAH-66  Comanche  Helicopter. 


The  Boeing  Defense  and  Space  Group's  Helicopter  Division  of  Philadelphia, 
Pennsylvania  and  United  Technologies’  Sikorsky  Aircraft  Division  of  Stratford, 
Connecticut  were  awarded  the  demonstration/validation  (DEM/Val)  phase  of  the  contract 
for  the  Comanche.  These  contractors  divided  the  aircraft  into  two  sections  for  design  and 
fabrication.  Boeing  has  responsibility  for  the  tail  section  of  the  aircraft.  Sikorsky  has 
responsibility  for  the  forward  portion  of  the  aircraft  fuselage,  to  include  the  integration  of 
both  sections. 

Currently,  the  first  prototype  is  undergoing  extensive  flight  testing  at  the  Sikorsky 
test  flight  facility  in  West  Palm  Beach,  Florida.  Focused  on  "envelope  expansion",  the 
aircraft  had  flown  approximately  91  hours,  reached  171  knots  in  forward,  level  flight, 
achieved  bank  angles  of  45°,  and  sideward  flight  of  70  knots  as  ofMay,  1998.  Prototype 
#2  is  scheduled  to  begin  flight  testing  in  1999. 
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B.  PROBLEM  DEFINITION 


During  flight  testing  the  empennage  of  the  aircraft  experienced  buffeting  where  the 
horizontal  and  vertical  tail  vibrated  excessively  in  certain  flight  regimes.  Observable  from 
both  the  chase  aircraft  and  on  flight  data  records  for  strain  gauges  and  accelerometers, 
these  high  buffet  loads  created  higher  than  anticipated  fitting  loads,  particularly  in  the  tail 
region,  and  resulted  in  vibrations  in  the  crew  stations  and  at  the  nose  cone  where  the 
targeting  sensors  are  located.  Significant  effort  has  been  devoted  to  identifying  the 
sources  of  excitation  and  the  nature  of  the  structural  response.  One  such  method  is  shown 
in  Figure  1 .2.  Here,  strips  of  material  were  placed  on  the  empennage  to  allow  for  better 
visualization  of  both  the  airflow  and  vibration  of  the  structure. 


Figure  1.2:  Comanche  Prototype  Test  Flight. 
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Boeing  and  Sikorsky  have  also  conducted  several  wind  tunnel  experiments  to 
better  quantify  the  flow  coming  from  the  aircraft  fuselage.  The  first  was  completed  in 
April  1997  at  NASA  Langley.  Two  follow-on  tests  were  executed  at  the  United 
Technologies’  Research  Center  in  Hartford,  Connecticut  and  then  Boeing  Defense  and 
Space  Group’s  Helicopter  Division  in  Philadelphia,  Pennsylvania.  These  tests  consisted  of 
changing  aerodynamic  surfaces  to  modify  the  flow  impinging  on  the  tail  from  the  fuselage. 
Additionally,  the  flight  test  aircraft  had  several  structural  modifications  to  stiffen  the  tail 
section;  notably  a  bracket  to  secure  the  horizontal  and  vertical  tail,  use  of  a  rigid  flap  on 
the  right  trailing  edge  of  the  horizontal  stabilizer,  and  the  adding  of  tip  weights  to  the 
horizontal  stabilizer. 

This  research  determines  the  location  and  magnitude  of  empennage  loads  caused 
by  the  random  excitation.  Because  the  nature  of  the  excitation  is  assumed  to  be  a  random 
process,  spectral  analysis  is  used.  To  obtain  the  loads,  a  three  step  process  is  utilized. 
First,  from  aircraft  differential  pressure  transducers  and  accelerometers,  the  spectral 
content  of  the  excitation  and  response  is  determined.  Then,  using  a  NASTRAN  model 
modified  to  replicate  the  flight  test  aircraft,  frequency  response  functions  are  obtained 
between  selected  points  on  the  aircraft's  tail  and  the  accelerometers.  Finally,  using  this 
information,  estimates  of  the  airloads  are  calculated.  Having  provided  information  on  the 
nature  of  the  excitation,  structural  modifications  can  then  be  explored  to  reduce  the 
response  of  the  airframe  to  the  excitation.  Potential  methods  for  accomplishing  this  may 
be  based  on  shifting  airframe  natural  frequencies,  increasing  structural  damping,  or  adding 
vibration  control  devices. 
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II.  HELICOPTER  VIBRATIONS 


A.  GENERAL 

As  with  all  helicopters,  the  RAH-66  experiences  vibrations.  These  forces  can 
significantly  impact  crew  effectiveness,  imaging  and  targeting  performance  and  airframe 
and  equipment  reliability  and  maintainability.  With  the  use  of  complex  digital  avionics 
systems,  there  is  a  greater  emphasis  on  vibration  reduction.  As  such,  the  analysis  of 
dynamic  loads  and  modification  of  structures  to  create  a  low  vibration  environment  is 
essential.  The  first  step  in  reducing  vibrations  is  to  identify  the  source  of  excitation.  Once 
that  is  known,  the  aircraft  can  be  modified  so  that  the  structural  frequencies  are  offset 
from  the  driving  frequencies. 

Helicopter  vibrations  fall  into  three  categories.  First,  there  are  vibrations  due  to 
rotor  excitation,  which  are  at  frequencies  that  are  integral  multiples  of  the  rotor’s 
rotational  speed.  Second,  vibrations  due  to  random  aerodynamic  excitation  where  the 
frequency  being  observed  is  a  natural  frequency  of  the  structure  being  excited.  And 
finally,  self-excited  vibrations,  such  as  flutter  and  ground  resonance.  [Ref.  2] 

1 .  Vibrations  at  Integral  Multiples  of  Rotor  Speed 
The  primary  source  of  helicopter  fuselage  vibrations  is  the  rotor  system.  Here, 
harmonics  of  aerodynamic  loads  on  the  blade  give  rise  to  vibratory  response  of  the  blade. 
Since  the  blade  is  restrained  at  the  root,  the  blade  responses  result  in  root  shears,  which 
feed  from  the  rotor  head  into  the  fuselage  as  vibratory  shears  and  moments.  As  the  forces 
go  from  the  rotating  to  the  fixed  fuselage  system,  the  rotor  system  tends  to  act  as  a  filter. 
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For  an  N  blades  rotor  system,  the  frequencies  which  will  filter  through  are  those  at  N  and 
2N/rev.  Since  the  amplitudes  of  the  lower  harmonics  of  blade  loading  are  greater  than  the 
higher,  the  N/rev  vibration  of  the  fuselage  is  most  critical.  It  can  be  shown  that  N/rev 
fuselage  vibrations  in  the  fixed  system  are  the  result  of  the  N-l,  N,  and  N+l/rev  vibratory 
response  of  the  blades  in  the  rotating  system. 

In  addition  to  the  N/rev  or  NP  vibration,  there  are  secondary  excitation  frequencies 
which  also  must  be  taken  into  account.  Any  unbalance  of  the  rotor  system,  such  as  blades 
out  of  track,  will  give  rise  to  1/rev  or  IP  excitation  of  the  fuselage.  Dissymmetry  of  blade 
dampers  may  cause  troublesome  2/rev  excitation.  Also  important  is  the  1  x  tail  rotor 
rotational  speed. 

Forward  speed  plays  a  critical  role  in  vibrations.  Figure  2.1  is  a  typical  plot  of  the 
variation  in  vibration  levels  as  a  function  of  airspeed.  The  increase  in  vibrations  levels  in 
the  30  to  45  knot  range  can  be  attributed  to  blade  vortex  interaction  due  to  the  proximity 
of  the  rotor  to  its  wake  and  the  resulting  non-uniform  induced  velocity  field. 


Figure  2. 1 :  Typical  Helicopter  Vibration  Levels  Versus  Airspeed.  After  Ref  [2], 
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In  the  high  speed  flight  regime,  the  increase  in  vibration  with  increase  in  airspeed  is 
attributed  to  two  sources.  First,  the  aerodynamic  excitation  due  to  impingement  of  the 
rotor’s  wake  on  the  helicopter’s  horizontal  empennage.  Second,  amplification  of  blade 
section  aerodynamic  loads  associated  with  increased  cyclic  pitch  accompanied  by 


increasing  differences  in  dynamic  pressure  per  blade  revolution.  Figure  2.2  shows  the 
non-distorted  rotor  wake  geometry  for  both  low  and  high  forward  flight  speeds. 


2.  Random  Aerodynamic  Excitation 

This  type  of  vibration  is  often  associated  with  the  turbulent  downwash  from  the 
rotor,  which  impinges  on  the  fuselage  and  its  components.  Similar  to  white  noise,  this 
excitation  contains  many  harmonics  that  can  excite  fundamental  fuselage  or  empennage 
modes.  Often,  an  interesting  characteristic  of  this  type  of  excitation  is  that  it  is  present 
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only  within  a  certain  range  of  airspeeds.  Successful  solutions  of  this  problem  have 
included  modifying  the  structure  excited  and  by  altering  the  flow  from  the  rotor. 

3.  Self-Excited  Vibrations 

A  self-excited  vibration  is  one  that  results  in  divergent  oscillations  where  the 
system  damping  is  negative.  In  the  case  of  negative  damping,  the  damping  force,  which  is 
now  a  driving  force,  does  positive  work  on  the  system.  The  work  done  by  this  force  is 
converted  into  the  additional  kinetic  energy  of  the  increased  vibration.  A  self-excited 
vibration  cannot  exist  without  an  external  source  of  energy.  In  the  case  of  flutter,  the 
airflow  provides  the  source  of  energy. 

B.  VIBRATION  SPECIFICATIONS 

Vibration  specifications  for  Army  aircraft  are  governed  by  Aeronautical  Design 
Standard  (ADS)  -  27  and  MIL-STD-810D  [Ref.  1],  Human  factors  vibration 
requirements  are  expressed  by  ADS-27  in  terms  of  an  intrusion  index  which  is  the  square 
root  of  a  weighted  sum  of  squares  of  four  vibration  components  in  each  direction.  The 
purpose  of  this  is  to  bring  together  the  twelve  most  significant  crew  vibration  components 
into  a  single  number  for  evaluation.  A  separate  specification  exists  to  capture  the  lower 
frequency  one  per  main  rotor  revolution,  or  IP,  vibration. 

The  changes  in  vibrations  that  result  for  airspeed  and  maneuvering  are  addressed  in 
ADS-27  by  defining  limits  in  terms  of  flight  regions.  There  are  three  defined  regions,  with 
the  first  being  normal  flight  up  to  cruise  speed,  then  maneuvers  having  a  duration  less  than 
3  seconds,  and  lastly  maneuvers  having  a  duration  greater  than  3  seconds.  There  is  also  a 
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region  for  gunfire  vibration.  The  specifications  for  the  RAH-66  crew  seat  are  listed  in  the 
Performance  Weapon  System  Specification  [Ref.  1]  and  are  shown  in  Table  2. 1 . 
Additionally,  Figures  2.3  and  2.4  display  the  currently  measured  intrusion  index  in  the 
Pilot’s  and  Copilot’s  crew  station  versus  the  specification. 


Intrusion  Index 

MR  IP  (in/sec) 

Flight  Region  I 

1.2 

.15 

Flight  Region  II 

3.0 

.30 

Flight  Region  III 

4.0 

.40 

Weapons  Firing  Increment 

+2.0 

N/A 

Table  2.1:  Crew  Seat  Intrusion  Index  and  MR1P  Vibration  Limits. 


Because  of  the  frequency  bands  that  are  used  to  calculate  the  intrusion  index,  any 
random,  broad  band,  excitation  occurring  between  main  rotor  IP  and  NP  is  not  totally 
captured.  Therefore,  the  intrusion  index  does  not  indicate  the  contribution  of  the  broad 
band  excitation  to  total  response. 
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Figure  2.4:  Intrusion  Index  -  Copilot’s  Station. 
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C.  COMANCHE  VIBRATION  SOURCES 

1.  Main  Rotor  Hub  Forces  and  Moments 

As  noted  earlier,  the  primary  source  of  helicopter  fuselage  vibrations  is  the  rotor 
system.  As  the  forces  and  moments  are  transmitted  from  the  rotating  to  the  fixed  fuselage 
system,  the  rotor  system  tends  to  act  as  a  filter.  It  can  be  shown  for  an  “N”  bladed  rotor 
system  that,  assuming  all  blades  see  the  same  loading  at  the  same  azimuth  position,  only 
NQ  frequency  will  be  seen  in  the  fixed  system.  This  NQ  frequency  in  the  fixed  system  can 
arise  from  either  N-l,  N,  or  N+l  excitation  in  the  rotating  system.  Figure  2.5  is  a 
schematic  of  the  hub  showing  the  mechanism  of  blade-fuselage  coupling. 


11 


With  the  coordinate  system  given  in  Figure  2.5,  forces  in  the  rotating  blade 
coordinate  system  can  be  related  to  forces  in  the  fixed  fuselage  using  the  following 
transformation: 
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where  [A,]  is  the  transformation  relating  forces  in  the  rotating  r|,  C,  system  to  forces  in  the 
fixed  X,  Y,  Z  system. 

Thus  for  the  five-bladed  Comanche  helicopter: 

•  4/rev  and  6/rev  flapwise  blade  root  shears  result  in  5/rev  hub  pitching  and 
rolling  moments  in  the  airframe. 

•  5/rev  flapwise  blade  root  shears  feed  into  the  fuselage  directly  as  5/rev  vertical 
forces. 

•  4/rev  and  6/rev  chordwise  root  shears  produce  5/rev  airframe  hub  forces  in  the 
fore  and  aft  and  lateral  directions. 

•  5/rev  chordwise  root  shears  results  in  5/rev  hub  yawing  moments. 
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For  the  RAH-66,  the  nominal  rotor  speed  is  354.9  rpm.  This  places  the  1/rev  or 
IP  excitation  at  5.915  Hz,  with  the  5P  blade  passage  frequency  occurring  at  29.575  Hz. 
[Ref.  3] 

2.  Aerodynamic  Excitation 

Aerodynamic  excitation  to  the  airframe  and  empennage  arises  from  several 
sources: 

•  Main  rotor  wake  impingement. 

•  Main  and  tail  rotor  downwash  (induced  flow)  impingement. 

•  Flow  from  around  the  main  rotor  hub  and  pylon  assembly  structure 

•  Flow  from  attached  aerodynamic  surfaces. 

•  Flow  from  protruding  irregular  fuselage  structure,  ranging  from  FABS 
(forward  avionics  bays)  to  as  small  as  a  canopy  or  door  handles. 

The  interaction  of  this  flow  is  a  complex  subject,  very  difficult  to  predict,  and  the 
subject  of  numerous  ongoing  research  projects.  Rotorcraft  can  have  problems  where 
rotor  airloads  affect  airframe  loads  and,  in  turn,  the  airframe  effects  rotor  airloads  and 
blade  dynamic  response.  Some  examples:  The  Sikorsky  H-3  had  a  lateral  tail  shake  that 
had  its  onset  at  about  120  knots.  The  problem  was  solved  by  installing  a  fairing  or 
“beanie”  on  top  of  the  main  rotor  hub  to  smooth  out  the  flow  from  this  region  which  was 
impinging  on  the  aft  fuselage. 

In  development  of  the  Boeing  YUH-61  and  Sikorsky  YUH-60  helicopters,  both 
companies  obtained  a  significant  improvement  in  N/rev  vibration  levels  by  “raising  the 
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mast”,  that  is,  increasing  the  separation  between  the  main  rotor  hub  and  the  fuselage.  In 
development  of  the  Hughes  YAH-64,  the  initial  T-tail  design  was  scrapped  and  a  movable 
horizontal  stabilator  installed  to  eliminate  an  excessive  vibration  problem  associated  with 
the  T-tail.  A  visual  summation  of  the  typical  sources  of  tail  excitation  is  shown  in  Figure 
2.6. 


Main  Rotor  Vortices 


Random  Excitation 
from  Flow  Separation 


Figure  2.6:  Typical  Sources  of  Tail  Excitation 

An  experimental  investigation  of  the  interactions  between  a  rotor  and  a  T-tail 
empennage  was  conducted  by  Lieshman  and  Moedersheim  [Ref.  4]  at  the  University  of 
Maryland  Rotorcraft  Center  in  1997.  Using  a  1/6  scale  four-bladed  helicopter  rotor  with 
generic  body  shape,  a  T-tail  empennage  was  placed  down  stream  of  the  rotor  hub  as 
illustrated  in  Figure  2.7.  Both  the  vertical  and  horizontal  tail  surfaces  used  the  NACA 
0012  airfoil  section  with  a  constant  chord  of  .203  meters.  The  horizontal  tail  was  located 
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in  two  different  positions,  a  high  and  low  plane  configuration,  with  zero  degrees  angle  of 
attack.  Additional  parameters  of  the  wind  tunnel  test  are  contained  in  Reference  4. 


Figure  2.7:  Schematic  of  Rotor/Body/Empennage  System.  From  [Ref.  4], 


Through  their  analysis,  they  drew  the  following  conclusions: 


•  Flow  behind  the  rotor  was  highly  asymmetric.  The  two  low  pressure  regions 
trailing  from  the  edges  of  the  rotor  disk  provided  evidence  of  the  wake  rolling- 
up  into  two  larger  vortex  bundles. 

•  At  low  advance  ratios  the  flow  over  the  horizontal  tail  was  stalled.  The  flow 
gradually  reattached  as  the  advance  ratio  was  increased.  Airloading  was 
asymmetric  primarily  due  to  the  difference  in  strength  of  the  two  trailing  vortex 
bundles  trailing  from  the  sides  of  the  rotor  disk.  In  some  test  conditions,  a 
sharp  trailing  edge  peak  pressure  was  observed  at  the  trailing  edge  of  the 
horizontal  tail.  This  indicated  the  potential  of  a  junction  vortex  originating  at 
the  junction  of  the  vertical  and  horizontal  tail. 

•  Measured  time-varying  pressures  showed  the  flowfield  near  the  horizontal  tail 
was  highly  unsteady.  Below  advanced  ratios  of .  1 5  the  unsteady  pressure 
response  was  low.  As  the  wake  transitioned  on  to  the  tail,  unsteady  pressures 
increased  significantly. 

•  Unsteady  pressure  responses  are  characterized  by  two  events,  unsteady  lift 
produced  by  the  rotor  wake  and  disturbances  caused  by  the  convection  of 
individual  vortices. 
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For  the  RAH-66,  the  rotor  wake  transitions  above  the  T-tail  between  40  and  60 
knots,  as  shown  in  Figure  2.5  [Ref.  3],  Beyond  this  airspeed  range  the  actual  wake  is  not 
impinging  on  the  horizontal  stabilizer.  From  Leishman’s  research;  however,  there  can  be 
convected  vortices  associated  with  the  passage  of  blades  that  can  induce  disturbances  and 
oscillating  loading  conditions. 


Figure  2.8:  Rotor  Wake  Distance  from  Horizontal  Stabilizer.  After  Ref.  [3], 

In  addition  to  the  rotor  vortices,  the  wake  behind  the  rotor  hub  constitutes  a  low- 
pressure  sink  that  can  draw  the  flow  off  of  the  upper  fuselage  and  produce  separation 
[Ref.  5],  The  frequency  content  of  this  excitation  is  often  determined  by  modeling  the 
rotor  hub  as  a  circular  cylinder  where  vortices  are  shed  as  a  Karman  vortex  street. 
Through  the  Strouhal  number,  a  non-dimensional  constant,  vortex  shedding  frequency  can 
be  related  to  freestream  velocity  and  the  diameter  of  a  two  dimensional  cylinder.  [Ref.  5] 
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The  relationship  is  given  by: 


where  fs  =  Vortex  shedding  frequency  (Hz) 

S  =  Strouhal  number 
U-  Freestream  velocity  (ft/sec) 

D  =  Cylinder  diameter  (ft) 

Determination  of  the  Strouhal  number  from  Schlichting  [Ref.  5:p:32]  requires 
interpolating  from  a  small  plot.  Using  the  results  from  a  wind  tunnel  experiment  done  at 
the  Naval  Postgraduate  School  on  a  two-dimensional  cylinder  by  Small,  Hebbar,  and 
Platzer  [Ref.  6],  the  Strouhal  number  and  vortex  shedding  frequency  can  more  easily  and 
accurately  be  determined.  In  this  experiment  [Ref.  6],  the  vortex  shedding  frequency  and 
Strouhal  number  were  calculated  for  Reynolds  numbers  ranging  from  2.7xl04  to  1.9xl05. 
From  this,  the  data  presented  in  Figure  2.9  was  generated.  The  Strouhal  number 
represents  the  slope  of  the  line  and  can  be  used  to  calculate  the  frequency  of  shedding 
from  a  know  diameter  and  given  velocity. 
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Figure  2.9:  Shedding  frequency  for  the  Quasi-two-dimensional  Cylinder.  After  Ref.  [6], 


For  the  Comanche,  at  100  knots  (168.78  ft/sec)  and  with  a  rotor  hub  pylon  of 
approximately  .75  meters  (2.46  ft),  the  calculated  vortex  shedding  frequency  is  1 1.827  Hz. 
As  the  airspeed  and  point  of  separation  changes,  the  shedding  frequency  will  be  effected, 
as  shown  by  the  equation  on  the  previous  page.  The  area  around  the  main  rotor  pylon  is 
actually  a  complex  array  of  surfaces  needed  to  reduce  radar  signature.  Separated  flow 
could  occur  from  any  of  these  surfaces,  with  varying  frequency  content  of  the  vortices.  In 
this  case,  more  broadband  vortex  shedding  will  occur. 

From  these  factors,  RAH-66  response  data  should  show  energy  associated  with 
main  rotor  5P,  29.75  Hz  and  if  a  rotor  track  and  balance  problem,  main  rotor  IP  at  5.95 
Hz.  There  could  also  be  energy  associated  with  separated  flow  from  the  main  rotor  pylon, 
centered  around  11-12  Hz,  and  some  broad  band  excitation  from  the  convected  rotor 
vortices  impinging  on  the  horizontal  tail. 
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III.  SPECTRAL  ANALYSIS 


Figure  3. 1  shows  the  time  history  of  lateral  excitation  at  the  pilot’s  heels.  When 
analyzed  over  the  foil  record  length  of  60  seconds,  the  data  is  classified  as  random.  Such 
random  time  functions  as  this  data  set  commonly  exist  in  nature;  a  typical  example  being 
pressure  gusts  aircraft  encountered  in  flight.  The  method  used  to  characterize  these 
functions  is  known  as  spectral  analysis,  and  a  discussion  of  the  theory,  from  Thomson 
[Ref.  7],  and  techniques  for  digital  signal  analysis  follows. 


Figure  3.1:  Time  History  of  Lateral  Acceleration. 


A.  TIME  AVERAGING 

Despite  the  irregular  nature  of  the  time  history  in  Figure  3.1,  most  “random” 
functions  exhibit  some  degree  of  statistical  regularity  and  by  applying  averaging  processes, 
we  can  characterize  the  phenomena.  In  observing  the  data,  each  record  is  referred  to  as  a 
sample  and  the  collection  of  samples  is  the  ensemble.  If  the  ensemble  average  is  taken  at  a 
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specific  time,  tj ,  and  then  at  another  time,  ti  +  t,  and  they  remain  the  same,  the  process  is 
said  to  be  stationaiy.  Similarly,  if  the  ensemble  averages  are  replaced  by  time  averages,  or 
averaging  each  sample,  and  the  averages  from  different  samples  remains  the  same,  the 
random  process  is  ergodic.  Key  assumptions  necessary  for  our  analysis  is  that  the  random 
response  phenomena  exhibited  by  the  aircraft  is  stationary  and  ergodic.  As  such,  statistical 
properties  can  be  derived  from  single  time  histories  of  sufficient  length.  [Ref.  7] 

To  derive  these  properties  a  review  of  statistics  is  required.  From  the  time 
histories,  we  can  define  the  expected  value  (mean): 

4*W]=lim-f  }*(')<* 

T  — *  Q 

For  the  case  of  discrete  variables,  the  expected  value  is  given  by  the  following  equation: 

«->•  00  fl  1=] 

We  can  then  apply  the  averaging  process  to  the  mean  square  value,  also  referred  to  as  the 
square  of  the  expected  value. 

T 

=\\mUx2dt 

T  — »cc  1  q 
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Linking  this  back  to  basic  statistics,  with  the  variance  defined  as: 


a2  =x2-(x)2 

°2  =\im\Vx-x)2dt 

T  — »CO  ■*  Q 

the  mean  square  value  (MSV)  is  simply  the  variance  when  the  data  has  an  expected  value 
(mean)  of  zero.  If  thought  of  in  terms  of  a  sin  wave,  or  multiple  harmonic  function,  where 
the  mean  is  removed  (detrended)  before  analysis,  then  the  variance  and  MSV  are  the  same. 

Knowing  these  statistical  properties  of  a  given  signal,  the  correlation,  or  similarity, 
between  different  signals  can  be  determined.  This  is  an  important  function  in  time  domain 
analysis  as  it  can  help  identify  signals  buried  in  noise  or  propagated  in  time  through  a 
structure.  Additionally,  it  provides  a  transition  into  Fourier  analysis  and  the  frequency 
domain,  where  a  structure’s  spectral  characters  can  better  be  analyzed.  The  correlation 
function  is  computed  by  multiplying  the  ordinates  of  two  records  at  each  time  and 
determining  the  average  value  by  dividing  the  sum  by  the  number  of  products.  The 
product  is  a  maximum  when  the  two  signals  are  very  similar  or  identical.  For  dissimilar 
records,  the  products  will  be  negative  and  positive,  so  the  sum  will  be  small.  If  two 
signals  are  identical  with  one  shifted  by  some  time,  we  can  compute  the  correlation  by 
taking  just  a  single  time  record  and  multiplying  it  by  itself,  with  the  specified  time  delay  i. 
This  is  called  the  autocorrelation  function,  and  is  defined  by  the  following  formula. 

1  % 

R(r)  =  E[x(t)x(t  +  t)]  =  (x(t)x(t  +  T)\  =  limv  f  *(0 +  T)dt 

r-> *  I  Jr/ 
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If  the  time  delay  t  =  0,  the  equation  reduces  to  the  mean  square  value  or  variance. 

R(  0)  =  x*  =  cr2 

If  the  two  records  are  truly  dissimilar,  we  call  their  product  the  cross  correlation.  It  is 
define  it  by: 

R„  (T)  =  E[x(t)y(t  +  t)]  =  (x(t)y(t  +  r)) 

1  % 

=  lim7  j  x(t)y(t  +  r)dt 

T-+ao  1  t  / 

~/2 

B.  FREQUENCY  DOMAIN  ANALYSIS 

Having  defined  averaging  processes  for  the  time  domain,  the  random  function  can 
be  transitioned  to  the  frequency  domain  [Ref.  7],  This  will  assist  in  relating  the  response 
data  back  to  the  sources  of  excitation  discussed  earlier. 

The  time  signal  x(t)  can  be  expressed  in  terms  of  its  Fourier  series: 

*(')  =  t,c„e"-' 

-oo 

Because  the  record  is  periodic  it  can  be  represented  by  the  Fourier  series: 


*(/)=Re2C.«^'  = 

n= 1 


-£(C.e*«  +C>inaV) 
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where  C„  a  complex  function  and  C„*  is  its  complex  conjugate.  It  can  then  be  proven  [Ref. 
7]  that  the  mean  square  value  of  the  function  is  given  by 


*2=ik,c; 

n-\  * 

From  this,  it  can  be  seen  that  the  MS  V  can  be  computed  in  either  the  time  or  frequency 
domain,  depending  on  which  technique  provides  better  analysis  of  the  data. 

The  mean  square  value  (MSV)  represents  contributions  from  the  entire  frequency 
band  of  the  signal.  If  only  a  certain  frequency  range  is  of  interest,  the  MSV  for  that  range 
of  frequency,  Af,  is  the  power  spectrum  (PS).  The  power  spectrum,  P(frJ,  is  then  just  the 
product  of  the  Fourier  coefficients  from  the  frequency  range  of  interest,^,  with  the  mean 
square  value  being  the  sum  of  the  power  spectrum  across  the  entire  bandwidth.  [Ref.  7] 

p<j,)=\c„c:  7=f 

n= 1 

If  the  power  spectrum  is  divided  by  the  frequency  interval,  Af,  the  result,  defined 
as  the  discrete  power  spectral  density(PSD)  S(frJ,  gives  the  contribution  to  the  power 
spectrum  from  an  even  more  narrow  frequency  band.  This  can  be  extended  for  a 
continuous  spectrum  as  the  limiting  case  of  S(fn)  as  A/— >  0. 


S(f„)  = 


pj/J  _  cmcrm 

A/  2A/ 


The  mean  square  value  can  then  be  computed  from  the  PSD  by 


H=1 
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The  root  mean  square  (RMS),  synonymous  with  the  standard  deviation,  can  be  found  by 
taking  the  square  root  of  the  mean  square  value  (MSV). 


RMS  =  4^ 


As  we  are  dealing  with  non-periodic  functions  over  discrete  time  intervals,  use  of 
the  Fourier  integral,  a  limiting  case  of  the  Fourier  series,  is  required.  The  Fourier  Integral 
and  Fourier  transform  of  x(t)  are  defined  respectively  as 


co  oo 

*(<)  =  I  X(f)e™df  X<J)  =  J 

-oo  -oo 

Using  Parseval’s  theorem,  the  MSV  of  the  discrete  signal  in  the  time  domain  can 
be  equated  to  the  Fourier  transform  of  the  signal  multiplied  by  its  complex  conjugate.  If 
Xi(t)  and  x 2(t)  are  real  valued  time  functions  and  their  Fourier  transforms  are  given  by 
Xi(f)  and  X2(f)  then 


J*,(0*2(0<#=  jx,(j)xl(/)d/ 


If  X\(t)  =  x2(t)  =  x(t)  the  MSV  and  PSD  can  be  expressed  in  terms  of  the  Fourier  transform 
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As  with  the  Fourier  coefficients,  the  MSV  and  RMS  of  the  signal  can  be 
determined  from  the  Fourier  transforms.  Additionally,  the  autocorrelation  function, 
symmetric  about  the  origin,  can  be  related  to  the  spectral  density  function  via  the  Weiner- 
Khintchine  equation.  [Ref.  7] 

oo 

S(f  )  -  2 J  R(t)  cos  Irfrch 

0 

The  cross-spectral  density  (CSD),  defined  as  the  Fourier  transform  of  the  cross¬ 
correlation  function: 

-oo 

Can  also  be  determined  by  multiplying  the  Fourier  transforms  of  the  two  signals  of 
interest. 

■*,(/> =  limi*'</)W) 

r-».  oo  1 

Given  all  of  these  tools,  the  one  of  major  interest  in  analyzing  the  flight  test  data  is 
the  power  spectral  density  (PSD).  This  function  of  frequency  gives  an  indication  of  the 
contribution  to  the  total  response,  in  our  case  accelerations,  due  to  a  discrete  frequency 
bandwidth.  Additionally,  the  cross-spectral  density  recaptures  the  phasing  information 
lost  when  the  Fourier  Integral  is  multiplied  by  its  complex  conjugate. 
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C.  DIGITAL  SIGNAL  PROCESSING 


Given  this  theory,  digital  signal  processing  methods  have  been  devised  that 
compute  the  PSD  from  discrete  time  histories.  Vectors  represent  the  time  histories.  In 
our  case,  the  sampling  rate  is  300  samples  per  second  and  record  length  is  60  seconds. 
Therefore,  the  vector  of  the  discrete  signal  has  18,000  data  points. 

To  convert  this  time  signal  into  either  a  PSD  or  CSD,  the  first  step  is  to  decide 
how  many  points,  N,  are  needed  in  the  Fast  Fourier  Transform  (FFT).  For  computational 
speed  and  accuracy,  FFTs  are  done  in  power  of  two.  The  number  of  points  also 
determines  what  the  effective  bandwidth  of  the  processed  information  will  be.  Bandwidth, 
Af,  is  defined  as  1/Time,  and  for  a  discrete  signal,  Time  =  N  x  At  (time  increment).  The 
time  increment.  At,  is  equal  to  1/sampling  rate,  or,  in  our  case ,  1/300  sec.  If  a  1024 
point  FFT  is  utilized,  Af =  300/1024  =  .293  Hz.  As  the  number  of  points  in  the  FFT 
increases,  the  bandwidth  decreases  and  you  get  finer  frequency  resolution. 

There  is  a  penalty  when  working  with  random  data  to  increasing  the  number  of 
points  in  the  FFT.  Based  on  N,  the  signal  is  divided  into  na  records  or  windows  and  has 
any  linear  trend  removed.  As  N  increases,  the  number  of  windows,  nj,  decreases  and  there 
is  less  averaging  occurring.  This  is  demonstrated  when  computing  the  expected  value. 
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To  suppress  leakage  of  the  signal  from  one  frequency  range  to  the  next,  time  windows  are 
used  that  taper  the  time  history  data.  This  eliminates  discontinuities  at  the  beginning  and 
end  of  the  record  to  be  analyzed.  A  commonly  used  window,  that  is  used  exclusively 
during  this  analysis,  is  the  Hanning  window.  It  is  given  by 

w(t)  =  .5  *  (1  -  cos(^-)) 

and  multiplied  by  the  N  data  points  of  each  nj  segment.  Use  of  the  window  does  reduce 
the  amplitude  of  the  processed  signal.  In  the  case  of  the  Hanning  window,  the  output  is 
.375  times  smaller.  A  correction  is  made  in  the  normalizing  of  the  PSD  to  account  for 
this. 

Variance  is  introduced  by  tapering  the  signal.  To  reduce  this,  the  n^  segments  are 
overlapped.  In  terms  of  the  Hanning  window,  the  overlapping  brings  the  "peaks"  of  each 
cosine  window  closer  to  each  other.  This  process,  basically,  takes  the  points  that  were  at 
the  end  of  one  n<i  segment,  where  the  window  is  near  zero,  and  makes  it  the  "main  area" 
for  the  next.  The  FFTs  are  then  computed  for  each  tapered,  overlapped  segment. 

From  the  FFTs,  the  PSD  can  be  determined  using  the  following  formulas: 

±(\fft\\  2*  E[<\FFjf)]*  At 

e[(\fft\)  ]  —  “  (norniwindow))1 

The  expected  value  is  the  sum  of  the  FFTs  of  each  segment  multiplied  by  its 
complex  conjugate,  or  the  magnitude  of  the  FFT  squared.  Either  way,  it  produces  a  real 
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valued  function  that  is  then  divided  by  the  number  of  segments,  n<j.  As  discussed  earlier, 
and  seen  here,  increasing  N  reduces  the  number  of  FFT  averages  in  the  expected  value. 

The  PSD  is  then  determined  by  "normalizing"  the  expected  value.  This  first 
involves  multiplying  by  2  as  the  PSD  represents  the  one-sided  autospectral  density 
function.  Then,  the  signal  is  multiplied  by  the  time  increment,  At  (or  dividing  by  the 
sampling  rate,  from  the  definition  of  At)  and  divided  by  the  norm  of  the  window,  squared. 
The  reasons  for  this  come  from  the  definition  of  the  FFT  and  by  canceling  several  terms. 
Bendat  and  Piersol  [Ref.  8:p.371]  provide  the  exact  proof;  however,  their  notation 
requires  dividing  the  signal  by  N  and  then  adjusting  for  the  attenuation  caused  by  a 
tapering  window.  The  norm  of  the  window  (its  length  in  N  dimensions)  accounts  for  both 
elements.  For  no  windowing,  computationally  done  using  a  boxcar  window,  the  norm2  is 
the  number  of  points  in  the  FFT,  1024.  For  a  Hanning  window,  the  norm2  =  1024  x  .375 
=  384.  As  discussed  earlier,  this  corresponds  to  the  attenuation  causing  by  using  a 
Hanning  window.  Figure  3.2  and  3.3  provide  a  visual  illustration  of  the  procedure  used  to 
compute  the  autospectral  density  function  (PSD). 

For  the  cross  spectral  density,  the  same  steps  are  followed  for  each  time  history. 
Then  the  FFTs  from  each  signal,  x(t)  and  y(t)  are  multiplied  by  each  other  and  the 
expected  value  is  a  complex  function.  E[FFT e]  is  normalized  using  the  same  method. 
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The  time  history  is  divided  into  segments.  For  computational  speed  and  accuracy,  each  segment 
consists  of  2n  data  points.  The  data  is  detrended  to  eliminate  the  stationary  component. 


A  tapering  window,  such  as  a  Hanning,  can  be  used  to  reduce  the  spectral  leakage.  To 
counter  the  variance  caused  by  tapering,  overlapping  time  histories  are  used. 


Hanning  Window 
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Figure  3.2:  PSD  Computational  Steps  -  Time  Domain. 
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Figure  3.3:  PSD  Computation  -  Frequency  Domain. 
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The  PSD  or  CSD  output  vector  is  of  length  N/2  with  a  frequency  band  from 
0:Af: (sampling  rate/2).  As  an  example,  with  N  of  1024  points  and  a  sampling  rate  of  300 
Hz,  the  maximum  frequency,  commonly  referred  to  as  the  Nyquist  frequency,  is  150  Hz. 
The  frequency  vector  begins  at  0  and  goes  to  150  Hz  in  increments  of  300/1024  or  .293. 
For  each  point  of  the  frequency  vector,  the  PSD  vector  has  a  corresponding  value.  This 
value  represents  the  energy  associated  with  that  frequency  band.  So,  in  this  example,  the 
first  value  of  the  PSD  vector  represents  the  energy  present  within  the  0  to  .293  Hz  range 
Units  of  either  the  PSD  or  CSD  are  (input  units)2/Hz.  In  the  case  of  accelerometer 
output,  the  PSD  units  are  g2/Hz.  With  units  like  this,  it  is  often  difficult  to  grasp  an 
understanding  of  what  the  information  really  means.  The  RMS  is  much  easier  to  interpret, 
which  has  the  same  units  of  the  input,  in  this  case,  g’s.  To  compute  the  RMS,  the  PSD 
vector  is  multiplied  by  the  frequency  increment,  Af.  This  vector  is  then  summed  within  the 
frequency  range  of  interest.  The  resulting  scalar  is  the  MSV  and  the  RMS  is  the  square 
root.  If,  for  example,  the  RMS  between  5-25  Hz  is  needed,  the  PSD  vector  is  multiplied 
by  Af  =  300/1024.  Then,  only  the  indices  of  the  vector  that  represent  energy  between  5- 
25  Hz  are  summed,  points  17-86  of  the  512  total  PSD  values.  The  square  root  is  the 
RMS  and  represents  the  total  energy  or  response  between  5  and  25  Hz.  Comparing  back 
to  the  time  domain,  this  RMS  value  equates  to  the  standard  deviation. 
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IV.  FLIGHT  TEST  DATA  ANALYSIS 


Time  histories  of  accelerometers  and  pressure  transducers  from  flight  tests  were 
analyzed  using  both  1024  and  4096  point  FFT.  For  computational  speed,  1024  point 
FFTs  were  utilized  for  final  analysis.  Several  different  tapering  windows  were  also  applied 
with  no  noticeable  difference.  The  excitation  frequencies  are  clearly  evident  and  are  an 
order  of  magnitude  higher  than  the  energy  of  the  surrounding  frequency  bands.  Therefore, 
even  when  no  side  lobe  reduction  techniques  were  utilized,  the  frequency  content  of  the 
response  was  pronounced.  To  remain  consistent  with  the  methods  used  by  Sikorsky,  a 
Hanning  window  with  67%  overlap  was  used.  This  did  reduce  the  spectral  leakage  and 
variance  associated  with  windowing  and  produced  a  more  accurate  PSD  value. 

A.  ACCELEROMETERS 

Data  records  at  60,  80,  100,  120,  and  140  knots  were  obtained  for  the 
accelerometers  listed  in  Table  4.1.  The  sampling  rate  for  all  signals  was  300  samples  per 
second  with  a  60  second  record  length.  A  complete  listing  of  all  test  flight  parameter  is 
contained  in  Appendix  A. 

The  PSD  of  each  accelerometer  through  the  range  of  airspeeds  was  computed.  In 
looking  at  the  energy  across  the  entire  frequency  spectrum,  0-150  Hz,  there  was 
significant  energy  from  0-30  Hz,  and  none  of  magnitude  beyond  that  range.  Figure  4. 1 
through  4.3  display  the  PSD  for  nose  cone  lateral,  pilot  heel  lateral,  and  copilot  heel 
vertical.  The  programs  used  to  compute  the  PSDs  and  CSDs  is  in  Appendix  B. 
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Accelerometers 

Mnemonic 

Measurement 

LNOSE:D:300 

Lateral  Nose 

VNOSE:D:300 

Vertical  Nose 

LONOSE:D:300 

Longitudinal  Nose 

LPHEEL:D:300 

Lateral  Pilot  Heel 

VPHEEL:D:300 

Vertical  Pilot  Heel 

LPF:D:300 

Lateral  Pilot  Floor 

VPF:D:300 

Vertical  Pilot  Floor 

LOPF:D:300 

Longitudinal  Pilot  Floor 

LPO  VRHD  :D :  3  00 

Lateral  Pilot  Overhead 

VCOPHEEL:D:300 

Vertical  Copilot  Heel 

LCOPHEEL:D:300 

Longitudinal  Copilot  Heel 

VCOPF:D:300 

Vertical  Copilot  Floor 

LCPO  VRHD  :D :  3  00 

Lateral  Copilot  Overhead 

LVPYLNT  :D :  3  00 

Lateral  -  Vertical  Pylon  Top 

VST  ABTLT  :D :  3  00 

Vertical  -  Hor  Stab  Left  Tip 

LOSTABTL:D:300 

Longitudinal  -  Hor  Stab  Left  Tip 

VSTABTRT:D:300 

Vertical  -  Hor  Stab  Right  Tip 

LOSTABTR:D:300 

Longitudinal  -  Hor  Stab  Right  Tip 

Table  4.1:  Accelerometers. 


As  previously  noted,  primary  excitation  for  a  5-bladed  rotor  occurs  at  5P,  which  is 
29.6  Hz  for  the  Comanche.  As  expected  for  this  case,  5P  response  centered  at  29.6  Hz 
dominates,  especially  at  higher  airspeeds.  In  contrast,  the  amplitude  of  the  other  frequency 
bands  is  not  discemable.  To  provide  more  insight  at  response  due  to  other  than  5P, 
Figures  4.4  through  4. 18  displays  the  accelerometer  response  from  5-25  Hz.  Within  this 
frequency  band,  much  greater  detail  in  the  response  data  can  be  discerned.  Given  on  the 
next  17  pages  are  plots  that  indicate  the  frequencies  of  excitation  along  with  the  expected 
increase  with  airspeed. 
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1  -3  PSD  Copilot  Heel  -  Vertical  Response 


Airspeed  Frequency 


Figure  4. 10:  Copilot  Heel  -  Vertical  5-25  Hertz. 
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ZH/76 


PSD  Pilot  Overhead  -  Lateral  Response 


Airspeed  Frequency 

Figure  4.12:  Pilot  Overhead  -  Lateral  5-25  Hertz. 
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PSD  Copilot  Overhead  -  Lateral  Response 
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Figure  4. 13 :  Copilot  Overhead  -  Lateral  5-25  Hertz. 
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g  /Hz 


PSD  Horizontal  Stab  Left  Tip  -  Vertical  Response 


AirsPeed  Frequency 


Figure  4.15:  Horizontal  Stab  Left  Tip  -  Vertical  5-25  Hertz. 
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g  /Hz 


PSD  Horizontal  Stab  Left  Tip  -  Longitudinal  Response 


Airspeed  Frequency 

Figure  4. 17:  Horizontal  Stab  Left  Tip  -  Longitudinal  5-25  Hertz. 
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From  a  review  of  the  plots,  there  are  three  frequency  bands  that  continuously 
display  energy.  Energy  in  this  context  is  based  upon  measured  vibration  response  data  as 
either  acceleration  in  g’s  or  PSD  in  (g’s)2/Hz.  The  first  band  is  centered  around  5-8  Hz, 
the  second  from  10-13  Hz,  and  the  last  22-25  Hz.  From  the  PSD  data,  the  RMS  for  each 
airspeed  was  calculated.  Figures  4.19  through  4.22  show  the  computed  RMS  from  the 
nose  cone,  pilot’s  station,  and  copilot’s  station.  The  RMS  including  and  not  including  5P 
is  plotted  along  with  both  lateral  and  vertical  response.  The  5P  vibration  adds  a  significant 
amount  to  the  total  response  and  produces  a  sharp  rise  in  accelerometer  readings 
acceleration  at  airspeeds  above  100  knots.  The  total  response  for  the  frequency  band  of 
5-25  Hz  begins  to  level  off  at  100  knots,  the  slope  being  much  less  than  the  RMS  with  5P 
included.  This  excitation,  not  including  5P,  is  less  prevalent  at  lower  airspeeds  and 
significantly  increases  from  80-100  knots.  Past  100  knots,  this  frequency  band  energy 
begins  to  level  off  while  5P  vibration  levels  increase  significantly.  By  plotting  both  lateral 
and  vertical  RMS,  comparisons  as  to  the  more  dominant  response  can  be  made.  Lateral  to 
vertical  ratios  typically  range  of  the  order  from  1.25: 1  to  2.0: 1 .  In  all  cases  except  for  the 
pilot’s  floor,  lateral  response  dominates.  To  capture  the  maximum  excitation  outside  of 
the  5P  range,  and  minimize  5P  response,  a  more  detailed  analysis  is  presented  of  the  PSD 
at  100  knots. 

Figures  4.23  through  4.28  present  the  PSD  for  100  knots.  Again,  both  lateral  and 
vertical  responses  are  plotted  together.  As  was  seen  in  the  RMS  calculations,  the  lateral 
response  is  dominant.  Additionally,  the  three  frequency  regions  identified  earlier  are 
evident. 
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Acceleration  (g's)  Acceleration  (g's) 


Root  Mean  Square  Nose  Cone 


Figure  4. 19:  Root  Mean  Square  -  Nose  Cone. 


Root  Mean  Square  Pilot  Heels 


Figure  4.20:  Root  Mean  Square  -  Pilot  Heels. 
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Acceleration  (g's)  Acceleration  (g's) 


Root  Mean  Square  Pilot  Floor 


Figure  4.21:  Root  Mean  Square:  Pilot  Floor. 
Root  Mean  Square  Copilot  Station 


Figure  4.22:  Root  Mean  Square  -  Copilot  Station. 
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g2/Hz  g  /Hz 


Frequency  (Hertz) 

Figure  4.27:  Horizontal  Stab  Left  Tip  -  100  Knots. 


PSD  Horizontal  Stab  Rt  Tip  -  1 024  Point  FFT 


Frequency  (Hertz) 


Figure  4.28:  Horizontal  Stab  Right  Tip  -  100  Knots. 
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To  quantify  the  amount  of  response  due  to  these  three  frequency  bands,  or  bins, 
RMS  values  for  the  discrete  ranges  of  5-8  Hz,  10-13  Hz,  and  22-25  Hz  were  calculated. 
Each  RMS  value  was  summed  and  compared  to  the  total  RMS,  not  including  5P  response. 
The  data  is  presented  in  Table  4.2.  Although  these  three  frequency  bins  represent  only  9 
Hz,  or  45%  of  the  total  bandwidth,  they  constitute  from  50%  of  the  total  response  to  as 
much  as  90%.  Comparing  the  bins,  the  more  dominant  frequency  range  was  10-13  Hz, 
with  60%  of  the  accelerometers  having  their  peak  response  in  this  region. 


Root  M 

ean  Square  - 1 00  Knots  -  g's  1 

Frequency  Bins 

%  Total 

RMS 

Accelerometer 

5-8  Hz 

10-13  Hz 

22-25  Hz 

Nose  Cone  Vertical 

0.0214 

0.0288 

51.87 

Nose  Cone  Lateral 

0.0448 

■MBM 

0.0306 

mBi 

Pilot  Heel  Lateral 

0.0352 

M 

0.0159 

IMBIilMi 

Pilot  Heel  Vertical 

0.0161 

0.0096 

57.09  1 

Pilot  Floor  Lateral 

0.0238 

0.0191 

0.0058 

WMMMl 

Pilot  Floor  Vertical 

0.0095 

mm 

0.0050 

64.88 

CPG  Heel  Vertical 

0.0054 

0.0103 

82.59 

CPG  Floor  Vertical 

0.0030 

0.0052 

HIM 

65.97 

Pilot  Overhead  Lat 

0.0215 

0.0816 

0.0868 

88.90 

CPG  Overhead  Lat 

0.0061 

0.0229 

gjKSi 

Vert  Pylon  Top  Lat 

0.1669 

0.1024 

87.02 

HorStabLt  Tip  Vert 

0.2928 

0.5035 

■M 

84.60 

HorStab  Rt  Tip  Vert 

0.4284 

0.5327 

81.18 

HorStab  Lt  Tip  Long 

0.2325 

0.3360 

75.25 

HorStab  Rt  Tip  Long 

0.2354 

■iEilll 

0.0987 

86.69 

Table  4.2:  RMS  of  Frequency  Bins  -  100  Knots. 
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B.  FLIGHT  TEST  PRESSURE  MEASUREMENTS 


The  differential  pressure  transducer  data  records  that  correspond  to  the 
accelerometer  data  are  listed  in  Table  4.3.  The  sampling  rate  remained  300  samples  per 
second  with  a  60  second  record  length. 


Differential  Pressures 

Mnemonic 

Location 

VERTDP2 

Differential  Pressure,  7%  chord,  30% 
up  Vertical  Tail 

VERTDP4 

Differential  Pressure,  7%  chord,  60% 
up  Vertical  Tail 

LHORZDP3 

Differential  Pressure,  7%  chord,  60% 
Span  Left  Horizontal  Tail 

RHORZDP4 

Differential  Pressure,  7%  chord,  60% 

Span  Right  Horizontal  Tail 

Table  4.2 

1:  Differential  Pressure  Transducers. 

The  frequency  content  of  the  differential  pressures  was  determined  and  compared 
to  the  response  PSDs,  findings  of  Leishman  and  Moedersheim  [Ref.  4],  and  the 
approximated  shedding  frequency  of  the  rotor  hub  pylon.  The  analysis  also  provided 
insight  into  which  areas  of  the  structure  are  excited.  This  aided  in  focussing  the 
NASTRAN  analysis,  deriving  the  transfer  functions  between  the  random  aerodynamic 
loading  and  the  response  locations,  and  correlation  with  mode  shapes.  This  will  be 
addressed  in  Chapter  VI.  Figure  4.29  shows  the  variation  of  differential  pressure  from  the 
left  horizontal  stabilizer  as  a  function  of  airspeed  from  60  to  140  knots. 
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Figure  4.29:  Differential  Pressure  -  Left  Horizontal  Stabilizer. 

As  with  the  response  data,  the  magnitude  of  the  5P  differential  pressure  (29.85  Hz) 
is  an  order  of  magnitudes  greater  than  excitation  in  any  other  frequency  band.  To  gain  a 
better  appreciation  for  the  vibratory  pressures  at  other  than  5P,  the  transducer  output  is 
plotted  for  the  5-25  Hz  range  in  Figures  4.30  through  4.33. 

Given  are  three  dimensional  plots  of  airspeed,  frequency  and  magnitude  of  the 
PSD.  The  asymmetry  between  the  left  and  right  horizontal  tail  is  expected  due  to  the 
angle  of  the  vertical  tail.  The  right  horizontal  tail  exhibited  pressure  spikes  at  the 
previously  noted  discrete  frequency  bands  of  5-8  Hz,  10-13  Hz,  and  22-25  Hz. 
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Airspeed 


Frequency 


Figure  4.30:  Differential  Pressure  -  Right  Horizontal  Tail 
.  n-3  7%  Chord  60%  Span  Lt  Horizontal  Stab 
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Figure  4.31:  Differential  Pressure  -  Left  Horizontal  Tail 


Airspeed 


Frequency 


Unlike  the  discrete  excitation  of  the  right  horizontal  tail,  the  left  horizontal  tail  has 
broad  band  excitation  from  11-19  Hz.  Both  span  locations  on  the  vertical  tail  show 
symmetry  in  the  shape  of  broad  band  excitation,  similar  to  that  of  the  left  horizontal. 

There  is  a  decrease  in  magnitude,  on  the  order  of  12: 1,  as  you  move  up  the  vertical  tail, 
from  the  30%  to  60%  span  location. 

Figures  4.34  through  4.35  are  the  PSDs  of  differential  pressures  at  100  knots. 

Both  horizontal  and  vertical  transducers  are  contained  in  the  same  plot  for  better 
comparison.  For  the  horizontal  tail,  the  magnitude  of  the  pressures  seen  on  the  left  is 
much  great  than  that  of  the  right.  The  energy  is  peaked  within  the  1 1 .5  to  14  Hz  band, 
with  an  additional  peak  at  18  Hz.  The  right  displays  its  peak  pressure  at  6  Hz,  which 
corresponds  with  the  IP  of  the  aircraft. 

The  vertical  tail,  at  both  the  30%  and  60%  span  location,  have  the  same 
frequency  distribution  as  the  left  horizontal,  1 1 .5  to  14  Hz.  As  seen  earlier,  the  magnitude 
of  the  pressure  continues  to  diminish  as  you  move  up  the  vertical  tail.  Although  the 
pressure  distribution  between  the  vertical  tail  and  left  horizontal  is  similar,  the  pressures 
seen  at  30%  span  of  the  vertical  are  3.5  times  greater  than  that  of  the  left  horizontal.  The 
pressures  at  the  60%  span  of  the  vertical  and  the  left  horizontal  are  of  equal  magnitude. 
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Figure  4.34:  Differential  Pressure  Horizontal  Tail  - 100  Knots. 
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Figure  4.35:  Differential  Pressure  Vertical  Tail  - 100  Knots. 
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Linking  the  pressure  distributions  back  to  response  data,  RMS  values  were 
calculated  for  the  discrete  frequency  bins  used  earlier.  The  unit  for  the  RMS  of  the 
pressures  is  psi.  Table  4.4  presents  both  the  magnitudes  of  pressures  as  you  move  up  the 
vertical  tail  and  along  the  horizontal  as  well  as  the  dominant  frequency  band. 


Root  Mean  Square  - 1 00  Knots  -  psi  I 

Frequency  Bins 

Pressure  Transducer 

5-8  Hz 

10-13  HZ 

22-25  Hz 

7%  Ch  30%  Sp  Vert 

0.0263 

0.0503 

0.0268 

7%  Ch  60%  Sp  Vert 

0.0182 

0.0206 

0.0184 

7%  Ch  60  %  Sp  Lt  Hor 

0.0121 

0.0268 

0.0135 

7%  Ch  60%  Sp  Rt  Hor 

0.0066 

0.0057 

0.0048 

Table  4.4:  RMS  of  Differential  Pressures  -  100  Knots. 


Understanding  the  flow  over  the  empennage,  even  for  the  limited  amount  of 
transducers,  assisted  in  focussing  the  analysis  of  the  structural  response.  The  following 
chapter  presents  the  structural  analysis  and  correlation  between  aerodynamic  forces  acting 
on  the  surface  and  response  at  accelerometers. 
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V.  STRUCTURAL  DYNAMICS 


The  dynamics  of  the  multiple  degree  of  freedom  (MDOF)  Comanche  structure 
govern  the  transfer  of  random  loads  to  fuselage  response.  In  order  to  understand  the 
MDOF  system,  the  single  degree  of  freedom  case,  also  called  a  single  modal  oscillator, 
must  first  be  studied.  From  these  solution  methods,  the  dynamics  of  the  complex  structure 
can  then  be  quantified  using  a  systematic  modal  analysis  approach.  Of  particular  interest  is 
the  response  of  the  system  to  harmonic,  or  frequency  dependent  excitation.  The  quantity 
that  describes  this  is  called  the  frequency  response  function  (FRF)  or,  from  the  Laplace 
domain,  the  transfer  function  evaluated  at  (jco).  Working  from  a  vibrations  perspective, 
the  ERF  for  a  single  DOF  system  will  be  computed  and  then  applied  to  a  MDOF 
Comanche  structure. 

A.  SINGLE  DEGREE  OF  FREEDOM  (SDOF)  OSCILLATORS 

The  response  of  a  single  degree  of  freedom  system  to  a  forced  harmonic  excitation 
can  be  determined  from  the  equation  of  motion  that  describes  the  system.  The  SDOF 
oscillator  is  characterized  by 

mx + cx  +  kx  =  F  cos  (fit) 

where 

m  Total  (modal)  mass  participating  in  mode  of  vibration 

k  Total  (modal)stiffhess  participating  in  mode  of  vibration 

c  Total  (modal)  damping  participating  in  mode  of  vibration 
F  Magnitude  of  forcing  function 
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The  equation  can  then  be  divided  by  m  and  is  in  the  generalized  form, 

c  k  F  F 

x  + — x  + — x  =  — cos (Qt)  or  x  +  2£ct)  x  +  cq2x  = — cos (Qi) 

m  m  m  m 

The  total  solution  is  comprised  of  two  parts: 

*(0  =  *homo(0  +  *par«(0 

The  homogeneous  solution,  also  known  as  the  transient  or  free  response  is 

*homo  (0  =  ***  (Cj  cos (o)dt)  +  C2  s\n(6)j)) 
where  the  constants  C/  and  C2  are  determined  by  applying  initial  conditions.  The 
particular  solution  is  found  for  a  specified  forcing  function.  In  the  case  of  a  cosine  forcing 
function,  the  solution  is 

*part(0  =  ^COS(Q/-O) 

and  the  system  responds  at  the  frequency  of  excitation,  Q. 

To  solve  for  the  amplitude  and  phase  of  the  response,  A  and  q>,  complex  algebra  is 
used  and  the  general  forcing  function  becomes 

F(t)  =  Fe)Ql  =  F(cos(Q/)  +  _/sin(fl/))  where F  -  — 

m 

With  this  definition,  the  general  equation  for  the  forced  response  becomes 

x  +  2  £conx  +  co2nx  =  Fe}0i 

and  if  the  assumed  solution  for  the  particular  solution  is 

(0  =  A** 
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the  necessary  derivatives  can  be  taken  to  find  A  in  terms  of  known  quantities 
(-Q2  +  2£D  a>J  +  co1n)A  =  F 


dividing  by  the  natural  frequency  squared,  con2  and  grouping  the  real  and  imaginary  parts, 

A  =  Z(Q)*A  =  F 


r  a2) 

+  2  C—J 

LI 

<»„  \ 

with  the  force  amplitude  redefined  as 

F  F  = 

— —  =  —  =>F 
a>l  k 

The  terms  inside  the  brackets  represent  the  complex  frequency  dependent  constant, 
Z(Q),  known  as  the  impedence,  which  is  a  function  of  the  frequency  Q.  The  equation  can 
also  be  written  in  the  form 

Converting  to  polar  form  and  dividing  magnitudes  and  subtracting  angles,  a  solution  for  A 

lz^(n)+jz^(n)]A=F 


is  obtained 


A  = 


A/Zreal(Q)2+Zim^(Q)2 
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The  modal  parameters  defined  in  our  equation  of  motion  can  be  substituted  in  to  yield 


A=- 


-e  ;<t>  where  O  =  tan' 
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The  particular  solution  is  then  defined  as 

xpan  (/)  =  \A\em~^  =  |^]((cos(3Q/  -  O)  +  j  sin(Q/  -  <D)) 

As  the  forcing  function  was  a  cosine  function,  only  the  real  part  is  taken  and  the 
general  solution  is 

The  amplitude  of  the  steady  state  response  is  given  by 

where  F  =  — 
k 

If  the  amplitude  is  then  normalized  by  the  static  displacement  F/k,  the  frequency  response 
function  (FRF)  is  found 
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The  FRF  can  then  be  used  to  relate  any  forcing  function  at  a  given  frequency  to  the 
output.  If  the  forcing  function  is  defined  as  F(Q)  and  the  output,  say  displacement,  is 
Y(Q),  then  the  relationship  between  the  two  is  given  by 

Y(Ci)  =  H(Q)F(Q) 

B.  MODAL  ANALYSIS 

Knowing  how  to  solve  for  the  FRF  of  a  SDOF  system,  these  methods  can  be 
applied  to  the  MDOF  systems  used  to  represent  complex  structures.  In  order  to  solve  the 
coupled  set  of  equations  associated  with  MDOF  systems,  they  first  must  be  decoupled  into 
single  DOF  modal  oscillators.  Each  modal  oscillator  is  then  solved  in  the  procedure  just 
described  and  then  combined  to  form  the  total  solution.  This  method  is  called  mode 
superposition  and  is  consistent  with  the  matrix  methods  currently  used  for  most  structural 
analysis.  [Ref.  9] 

If  the  MDOF  structure  is  defined  as  a  set  of  NxN  coupled  equations  represented 
by 

[M]{ji}+[C]!i}+[KKl}={F(t)} 

The  solution  to  this  equation  is  first  found  by  solving  the  homogeneous  eigenvalue 
problem  with  no  damping,  C.  Substituting  an  assumed  solution  of  the  form  {x}={<( >}Qe'M 
and  dividing  by  the  mass,  the  eigenvalue  problem  is  defined  as 
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|m]"‘[K]-(02[i]|(o}={o} 


The  non-trivial  solution  of  the  eigenvalue  problem  produces  a  vector  of  length  N 
that  contains  the  natural  frequencies  squared,  con2  with  units  of  (rad/sec)2.  The  mode 
shapes,  or  eigenvectors  {<})},  can  be  easily  solved  for  and  combined  into  the  modal  matrix. 

The  modal  matrix,  although  not  an  orthogonal  matrix,  is  orthogonal  with  respect 
to  the  mass  and  stiffness  matrix.  As  such,  the  mass  and  stiffness  matrices  can  be 
transformed  by  the  modal  matrix  to  diagonal  matrices  that  decouple  the  equations  of 
motion  and  yield  N  -  SDOF  modal  oscillators.  The  coordinate  system  must  also  be 
transformed  to  modal  coordinates,  q 

W-MM 

Pre-multiplying  by  the  transpose  of  the  modal  matrix  and  substituting  for  x,  the  equation 
of  motion,  minus  the  damping  becomes 

Defining  the  modal  mass,  stiffness,  and  force  as 

[4.]t[mI«.]=[mJ;  [<*>]T[kI«>]=  [k];  {?}=  WiFi 
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The  modal  equation  of  motion  is 

[m|?}+[£]M={f} 

These  diagonal  matrices  now  represent  the  N  modal  oscillator  SDOF  equations  of  motion. 
A  solution  for  each  can  easily  be  obtained. 

Damping  is  re-introduced  on  a  modal  basis.  As  each  differential  equation  is 
solved,  the  damping  associated  with  that  particular  natural  frequency,  or  mode,  is  applied. 
This  is  consistent  with  the  techniques  used  in  parameter  identification  where  damping 
ratios  are  determined  on  a  modal  basis.  The  response  in  physical  coordinates  is  found  by 
multiply  each  modal  solution  by  the  modal  matrix 

H-MM 

If  interested  in  the  frequency  response  function  (FRF)  between  only  certain  points, 
i  and  j,  on  an  MDOF  system,  the  following  method  is  used 

"p"  modes  t&'-cb'- 

EL(Q)=  V  — - - 

rti  ^-Q2+j2Cra>rQ 

Once  the  FRF  is  known  for  those  points,  the  magnitude  of  the  response,  A,  due  to 
a  harmonic  excitation  of  amplitude  F  and  at  a  driving  frequency  of  Q,  can  be  determined 
by  multiplying  the  FRF  and  F. 
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C.  RIGID  BODY  MODES 


Unconstrained  systems  have  rigid  body  modes.  This  occurs  when  the  structure 
translates  or  rotates  without  any  deformation.  Numerically,  this  appears  in  the  solution  as 
natural  frequencies,  con,  of  zero.  The  reason  for  this  is  based  on  the  relationship  between 
the  stiffness  matrix,  K,  and  the  strain  energy  of  the  system.  If  the  MDOF  system  moves 
without  deforming,  no  strain  energy  is  produced.  The  stiffness  matrix,  K,  becomes 
positive  semi-definite  with  a  determinate  of  zero.  As  such,  the  solution  to  the  eignevalue 
problem  has  as  many  trivial  solutions,  zero,  as  there  are  rigid  body  modes.  For  an  aircraft, 
there  are  six  rigid  body  modes,  which  is  the  case  for  the  dynamics  of  any  body  in  three- 
dimensional  space. 

D.  FINITE  ELEMENT  MODELING 

To  determine  the  frequency  response  functions  using  modal  analysis,  the 

Comanche  structure  was  modeled  using  the  software  package  NASTRAN.  The  National 
Aeronautics  and  Space  Administration  (NASA)  initially  development  NASTRAN  in  the 
1960s.  NASTRAN  is  an  acronym  for  NASA  STRuctural  ANalysis.  Originally  written  in 
FORTRAN,  it  uses  the  finite  element  method  (FEM)  to  analyze  structural  models  [Ref. 
10],  The  foundation  of  this  is  the  discretization  of  the  structures  stiffness  and  mass.  The 
FEM  provides  the  basis  for  algorithms  that  can  efficiently  analyze  complex  structures  such 
as  the  Comanche.  A  detailed  discussion  of  the  Finite  Element  method  and  its  relationship 
to  structural  dynamics  is  contained  in  Reference  1 1 .  NASTRAN  version  70  was  used  for 
this  analysis. 
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The  results  generated  using  NASTRAN  were  viewed  using  PATRAN  7.0,  an 
integrated  computer-aided  engineering  (CAE)  tool.  The  menu-driven  graphical  interface 
makes  model  manipulation  easy  and  aids  in  visualization  of  mode  shapes  and  model 
geometry. 
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VI.  RAH -66  NASTRAN  ANALYSIS 

A  NASTRAN  finite  element  model  of  the  RAH-66  was  provided  by  Sikorsky.  This 
model,  yrahjl57_reinf_2kg.dat,  was  modified  by  engineers  at  Boeing  to  replicate  the 
configuration  of  the  prototype  aircraft.  The  name  of  the  model  is  very  important  as 
Sikorsky  manages  many  configurations  and  models,  each  with  different  stiffness  and  mass 
properties.  The  model  used  in  this  analysis  included  additional  stiffening  to  the  tail  of  the 
aircraft  and  2.2  kg  tip  weights  on  both  sides  of  the  horizontal  tail.  With  this,  the  natural 
frequencies,  mode  shapes,  and  frequency  response  functions  (FRF)  could  be  determined 
that  most  closely  represented  the  actual  aircraft.  This  was  critical  as  these  FRFs  were  later 
used  with  accelerometer  and  pressure  transducer  data  to  compute  estimates  of  the  loads. 

A  PATRAN  representation  of  the  NASTRAN  model  is  shown  in  Figure  6.1. 

A.  NATURAL  FREQUENCIES 

Hover  frequency  response  testing  of  the  test  flight  aircraft  was  conducted  at  West 
Palm  Beach  to  support  aeroservoelastic  analysis.  The  purpose  was  to  identify  actual 
natural  frequencies  and  modal  damping.  Through  flight  control  frequency  sweeps  from  5- 
1 5  Hz,  actual  fuselage  natural  frequencies  and  percent  critical  damping  were  determined 
for  selected  modes. 
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Figure  6. 1 :  RAH-66  NASTRAN  Model. 


Rather  than  adjusting  the  stiffness  and  mass  properties  of  the  NASTRAN  model  to 
return  the  actual  natural  frequencies,  engineers  at  Boeing  demonstrated  that  if  the  mode 
shapes  from  the  original  NASTRAN  model  were  used  along  with  adjusted  natural 
frequencies,  the  resulting  displacements  and  accelerations  matched  actual  test  flight 
response  data.  The  method  used  to  accomplish  that  was  to  first  run  the  NASTRAN 
analysis  deck.  Then,  during  a  restart  using  DMAP  commands,  the  eigenvalue  matrix 
(representing  NASTRAN  natural  frequencies)  was  altered  to  actual  frequencies.  The 
resulting  accelerations  from  application  of  harmonic  loads  closely  matched  actual  response 
due  to  harmonic  sweeps  of  the  tail  fan.  Table  6. 1  displays  the  NASTRAN  natural 
frequencies  and  actual  natural  frequencies  and  damping,  for  selected  modes,  as  determined 
from  flight  test. 


Aft  Fuselage  Lat/Torsion 

6.00 

Fuselage  Vertical  Bending 

7.70 

Pilot  Vert/Pitch 

Pilot  Lat/Y aw 

Horiz  Stab  Roll 

9.50 

Copilot  Vert/Pitch 

Copilot  Lat/Y  aw 

Horiz  Stab  Yaw 

11.60 

Tail  Rotor  Drive  Shaft  Mode 

2nd  Fuselage  Vertical 

Quill  Shaft 

Tail  Rotor  Drive  Shaft 


Fuselage  Head  Vertical 


2nd  Fuselage  Lateral 


Pilot  Seat  Lateral  Roll 


Main  Xmsn  Lat/Roll 


Tail  Rotor  Drive  Shaft 


Tail  Rotor  Drive  Shaft 


Ammo  Ba 


Copilot  Lat/Roll 


Main  Xmsn  Lat/Roll 


Fuselage  Mode 


Pilot  Seat 


Stab  Vertical  Bendin 


Tail  Rotor  Driveshaft 

Head  Lateral 

Tail  Rotor  Driveshaft 

Copilot  Sidearm 

Weapons  Door 

Fuselage 

Weapons  Door 
ECU  Mode 


24.74  1  Canopy _ 

Table  6. 1 :  Aircraft  Natural  Frequencies. 
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B.  MODE  SHAPES 


Selected  mode  shapes  for  free  response  from  the  original  NASTRAN  solution  are 
contained  in  Appendix  C.  Again,  as  the  flight  test  derived  natural  frequencies  were  only 
slightly  different  from  the  NASTRAN  frequencies,  it  is  assumed  that  the  NASTRAN 
mode  shapes  are  accurate.  Three  critical  mode  shapes  are  shown  in  Figure  6.2  through 
6.5.  The  first  two  mode  shapes  were  identified  in  test  flight  and  have  large  lateral 
response  in  frequency  bins  previously  discussed.  The  third  mode  shape  is  outside  of  the 
frequency  range;  however  it  correlates  well  with  the  response  seen  in  the  22-25  Hz  range. 

Figure  6.2  is  the  aft  fuselage  lateral  torsion  mode.  This  motion  is  consistent  with 
the  accelerometer  output  in  the  5-8  Hz  bandwidth.  Referring  to  Table  4.2,  correlation  can 
be  made  between  the  flight  test  output  and  the  free  response.  There  is  large  response  in 
this  frequency  bin  from  the  accelerometer  at  the  top  of  the  vertical  stabilizer.  This  same 
motion  occurs  in  the  6  Hz  mode  shape. 

The  horizontal  stabilizer  yaw  mode  at  1 1 .6  Hz,  displayed  in  Figure  6.3,  has 
significant  lateral  response  at  the  base  of  the  vertical  tail  and  at  both  the  pilot  station  and 
nose  cone.  The  RMS  of  accelerometer  output  from  both  the  pilot  station  and  nose  cone 
have  the  largest  magnitude  in  the  10-13  Hz  bandwidth. 

There  is  no  mode  shape  within  the  22-25  Hz  bandwidth  that  correlates  well  with 
the  flight  test  data.  Figure  6.4  shows  the  horizontal  stabilizer  bending  mode  at  26.02  Hz. 
In  the  baseline  NASTRAN  model,  yrahj  1 70,  the  frequency  is  22.08  Hz.  This  response  is 
consistent  with  flight  data  and  is  characterized  by  vertical  motion  of  the  horizontal 
stabilizer  and  lateral  response  of  the  CPG  station  and  pilot  overhead. 
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Figure  6.2:  Aft  Fuselage  Lateral  Torsion  Mode  Shape. 
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iw  Mode,  1 1.60  Hertz 


Figure  6.3:  Horizontal  Stabilizer  Yaw  Mode  Shape. 
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Figure  6.4: 


ibilizer  Vertical  Bending  Mode. 
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C.  FREQUENCY  RESPONSE  FUNCTIONS 


Frequency  response  functions  were  determined  between  the  horizontal  and  vertical 
tail  and  accelerometers.  To  compute  the  FRFs,  accelerometer  locations  had  to  be 
correlated  to  NASTRAN  grids.  Sikorsky  provided  a  list,  contained  in  Appendix  D,  of 
accelerometer  locations  and  suggested  NASTRAN  grids.  A  partial  list  of  the  NASTRAN 
grid  points  used  is  shown  in  Table  6.2. 


Accelerometer/; 

Correia? 

NASTRAN 

don 

Location 

NASTRAN  Grid 

Nose  Cone 

61641 

Pilot  Heel 

6505 

Pilot  Floor 

7246 

Copilot  Heel 

7608 

Copilot  Floor 

8948 

Pilot  Overhead 

7514 

Copilot  Overhead 

8614 

Vertical  Pylon  Top 

32341 

Hor  Stab  Left  Tip 

41108 

Hor  Stab  Right  Tip 

48108 

Table  6.2:  Accelerometer  Locations  in  NASTRAN  Model 


A  one  pound  force  was  applied  at  these  nodes  with  a  frequency  sweep  of  5-25  Hz. 
Accelerations,  in  g’s  and  in  all  three  directions,  were  found  at  each  accelerometer  location. 
The  FRFs  were  written  to  a  pch  file,  converted  into  matrices  compatible  with  MATLAB, 
and  stored  as  .mat  files  for  use  in  solving  for  PSDs  and  CSDs  of  loads. 
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1. 


Horizontal  Tail 


The  NASTRAN  nodes  along  the  quarter  chord  of  the  horizontal  tail  are  shown  in 
Figure  6.5.  Only  one  loading  condition,  loads  applied  in  the  z  direction,  was  used.  In 
addition  to  determining  the  FRF  at  each  grid  point  along  the  quarter  chord  and  the 
accelerometer  locations,  a  chordwise  investigation  was  conducted.  At  five  different  span 
locations,  the  FRF  was  found  for  each  node  along  the  chord.  This  provided  no  additional 
information  as  the  FRF  varied  minimally  from  leading  edge  to  trailing  edge.  Therefore,  it 
is  not  included. 


Figure  6.5:  Horizontal  Tail  Nodes  and  Loading  Condition. 
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Figures  6.6  and  6.7  show  the  FRF  between  the  horizontal  stabilizer  and  the  pilot 
heels  lateral  and  vertical.  Although  the  FRF  was  determined  from  5-25  Hz,  the  magnitude 
beyond  20  Hz  was  extremely  small.  For  better  resolution  in  the  plot,  the  axis  was  reduced 
to  5-20  Hz.  The  span  axis  represents  every  node  along  the  quarter  chord,  with  the  right 
tip  defined  as  0%  span  and  the  left  tip  representing  100%.  As  expected,  the  peak  FRFs 
are  aligned  with  the  natural  frequencies. 

In  Figure  6.6,  the  lateral  FRF  has  the  greatest  magnitude  at  6  Hz  on  the  right  tip  of 
the  stabilizer.  All  other  FRFs  are  less  than  this  value.  As  the  load  was  moved  along  the 
tail,  the  magnitude  continued  to  decrease  until  a  minimum  at  mid-span.  This  is  expected 
as  loads  directly  over  the  vertical  tail  induce  very  small  moments,  and  the  lateral  motion 
throughout  the  fuselage  reduces.  Continuing  out  the  left  horizontal,  the  magnitude 
continues  to  increase;  however,  never  reaches  the  amplitude  from  the  right  side. 

Figure  6.7  shows  the  FRF  between  the  applied  load  and  the  vertical  acceleration  at 
the  pilot  heels.  The  FRF  does  have  symmetry  across  the  tail;  however,  the  magnitude  is 
substantially  less  than  the  lateral. 

Figures  6.8  and  6.9  display  the  FRF  between  the  same  points  but  at  the  discrete 
frequencies  of  6  and  1 1 .6  Hz.  The  magnitude  was  normalized  to  the  maximum  amplitude 
FRF.  Comparisons  can  now  be  more  easily  drawn  for  relative  values  of  each  FRF.  These 
plots  indicate  the  significant  role  the  right  side  of  the  horizontal  tail  plays  in  transferring 
energy  at  both  6  and  11.6  Hz. 
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x  10 °  Horizontal  Tail  -  Pilot  Heel,  Vertical  Response 
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Figure  6.7:  FRF  Horizontal  Tail  to  Pilot  Heel  -  Vertical  Response. 
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%  Span  -  Right  Tip  to  Left  Tip 


Figure  6.9:  FRF  Horizontal  Tail  to  Pilot  Heel  -  Vertical. 


2. 


Vertical  Tail 


The  loading  condition  and  NASTRAN  nodes  along  the  quarter  chord  of  the 
vertical  tail  are  shown  in  Figure  6.10.  Again,  only  one  loading  condition,  loads  applied  in 
the  y  direction,  was  used.  As  with  the  horizontal  tail,  a  chordwise  investigation  was 
conducted.  Similarly,  the  FRF  chordwise  at  five  span  locations  varied  only  slightly. 

Those  results  are  not  presented. 


Figure  6. 10:  Vertical  Tail  Nodes  and  Loading  Condition. 
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Figures  6.11  and  6.12  show  the  FRF  between  the  vertical  stabilizer  and  the  pilot 
heels  lateral  and  vertical.  As  with  the  horizontal,  no  significant  response  level  was  evident 
in  the  20-25  Hz  range.  The  span  axis  represents  every  node  along  the  quarter  chord,  with 
the  base  of  the  tail  defined  as  0%  span  and  the  top  representing  100%. 

In  Figure  6. 1 1,  the  lateral  FRF  of  greatest  magnitude  is  at  6  Hz.  This  occurs  at  the 
top  of  the  vertical  tail.  All  FRFs  decreased  as  the  load  was  moved  down  the  tail,  except 
at  1 1 .6  Hz.  The  FRF  for  this  frequency  reached  a  maximum  at  the  base  of  the  vertical  tail. 

Figure  6. 12  shows  the  FRF  between  the  applied  load  and  vertical  acceleration  at 
the  pilot  heels.  It  is  consistent  with  the  information  from  the  previous  plot.  Magnitudes  at 
all  but  1 1.6  Hz  decrease  as  the  load  was  moved  down  the  tail.  The  amplitudes  are 
substantially  less  than  lateral. 

Figures  6.13  and  6.14  display  the  FRF  between  the  same  points  but  at  the  discrete 
frequencies  of  6  and  1 1 .6  Hz.  As  with  the  horizontal  stabilizer,  the  magnitude  was 
normalized  to  the  FRF  of  maximum  amplitude,  which  occurred  on  the  right  tip  of  the 
horizontal  stabilizer. 

The  top  of  the  vertical  has  the  same  magnitude  at  6  Hz  as  does  the  right  tip  of  the 
horizontal.  The  base  of  the  vertical  tail  has  the  dominant  FRF  at  1 1.6  Hz.  Comparing  this 
with  the  pressure  transducer  data,  the  right  horizontal  tail  had  the  most  energy 
concentrated  in  the  5-8  Hz  bandwidth.  This  is  where  the  6  Hz  FRF  is  a  maximum.  For 
the  base  of  the  vertical  tail,  the  pressure  energy  was  broadband  from  1 1-19  Hz.  The  1 1 .6 
Hz  FRF  is  a  maximum  at  this  location. 


91 


g 's/lb 


x  io'3  Vertical  Tail  -  Pilot  Heel,  Vertical  Response 


Span 


10  12.5 

Freauencv 


Figure  6. 12:  FRF  Vertical  Tail  to  Pilot  Heel  -  Vertical  Response. 
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Figure  6.13:  FRF  Vertical  Tail  to  Pilot  Heel  -  Lateral. 


Figure  6. 14:  FRF  Vertical  Tail  to  Pilot  Heel  -Vertical. 


VII.  RESULTS  AND  ANALYSIS 


A.  SOLUTION  METHOD 

Having  obtained  the  PSDs  of  the  response  and  the  structural  frequency  response 
functions,  a  solution  can  be  obtained  for  the  spectral  content  of  the  applied  loads.  The 
method  for  accomplishing  this  is  provided  by  Bendat  and  Piersol  [Ref.  8:p.  240],  A 
summary  of  this  method  follows. 

X  is  defined  as  a  column  vector  representing  the  Fourier  transform  of  q  input 
records,  and  Y  the  Fourier  transform  of  k  output  records; 


x  = 

*2 

Y  = 

?2 

/*_ 

X*,  Y*  =  complex  conjugate  of  the  vectors  X,  Y 
X',  Y'  =  transpose  of  the  vectors  X,  Y 

The  input  and  output  spectral  density  matrices  can  then  be  defined  as 


=  —  is{x*X'}=  input  spectral  density  matrix 


£{y*Y'}=  output  spectral  density  matrix 


95 


The  notation  used  by  Bendat  and  Piersol  for  and  G„  simply  means  that  the 
PSDs  and  CSDs  of  each  input  and  output  record  are  contained  with  these  matrices. 
Another  method  for  representing  these  terms,  and  useful  for  programming  is 


Both  Gxx  and  Gyy  are  Hermitian  matrices,  with  Gy  =  Gy,  for  all  i  and  j.  Each 
element  in  both  matrices  represents  either  the  PSD  or  CSD  for  the  inputs  and  responses. 
When  the  subscripts  are  the  same,  as  with  Gn,  the  value  is  a  real  number  and  represents 
the  PSD  for  the  load  at  point  one.  The  other  elements  of  the  matrix  are  complex,  such  as 
Gn  which  is  the  CSD  between  the  loads  at  point  one  and  two.  The  same  is  true  for  the 
response  matrix,  with  the  diagonal  being  the  PSDs  for  responses  and  the  off  diagonal 
terms  are  the  CSDs  between  responses. 

Finally,  the  FRF,  Hxy,  where  x  defines  the  input  and  y  is  the  output  is  given  by 


• 

■ 

- 
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Depending  on  the  number  of  inputs  and  outputs  used  in  the  solution  set,  the  FRF  matrix 
may  not  be  square.  Using  this  definition  of  Hxy,  it  follows  that  the  output  Y  is  related  to 
the  input  X  by  the  relationship 

Y  =  H^,X  where  is  the  transpose  matrix  to  H  ^ 

Having  defined  these  terms,  the  following  derivation  is  required  to  arrive  at  a 
solution  for  Gxx- 


y  =  h;x 

Y'  =  (H^X)'  =  X'Hjy 

Y*  =  (H^X)*=H^X* 

Y*Y  =  (H'iy*X*)(X'Hiy) 

Taking  the  expected  value  of  both  sides  and  multiplying  by  (2/T)  gives 

g„=(h;-)g=(hw) 

It  follows  that 

Assuming  the  required  inverse  matrices  exist,  a  solution  can  be  obtained  for  the 
spectral  content  of  the  applied  loads.  In  the  case  of  unequal  inputs  and  responses,  where 
Hxy  is  not  square,  the  solution,  Gxx,  is  found  using  gaussian  elimination.  Conditioning  of 
the  matrices  becomes  critical  during  the  computations. 
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B.  MATLAB®  PROGRAM 


Using  this  formulation,  a  MATLAB®  code  was  written,  contained  in  Appendix  E, 
that  solved  for  G^  First,  the  user  specified  the  frequency  range  of  interest.  This  was 
limited  from  5-25  Hz  because  of  the  range  of  frequency  sweeps  use  to  determine  the 
frequency  response  functions.  Then,  NASTRAN  grid  points  were  selected  for  application 
of  forces.  Finally,  response  points  to  be  used  in  the  solution  set  were  chosen. 

The  program  looks  up,  on  a  frequency  by  frequency  basis,  the  stored  FRFs 
between  force  application  points  and  response  locations.  This  created  the  Hxy  matrix. 

The  Gyy  matrix  was  found,  again  on  a  frequency  by  frequency  basis,  from  the  selected 
response  points.  PSDs  and  CSDs  for  each  set  of  accelerometer  flight  test  data  was  stored 
in  .mat  files.  This  allowed  retrieval  for  use  in  the  appropriate  Gyy  matrix.  Because  the 
frequency  interval  on  FRFs  was  .05  Hz  and  on  the  1024  point  FFT  PSDs  was  .293  Hz, 
cubic  spline  interpolation  was  used  to  match  the  data  sets. 

The  solution,  Gxx,  represented  the  PSDs  and  CSDs  of  the  applied  loads  at  discrete 
frequencies.  The  frequency  band  was  dependent  on  what  the  user  had  selected  for  the 
upper  and  lower  limit.  The  interval  between  frequencies,  or  effective  bandwidth  was 
based  on  Gyy.  Hence  for  this  output,  with  a  1 024  point  FFT,  the  interval  was  .293  Hz. 
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C.  RESULTS 


Numerous  cases,  with  different  locations  of  forces  and  reponses,  were  computed. 
Some  generalities  included  all  forces  placed  along  the  quarter  chord  of  either  the 
horizontal  or  vertical  tail.  These  points  are  shown  in  Figures  6.5  and  6. 10  and  listed  in 
Table  7.1.  Additionally,  the  only  response  data  used  was  from  a  100  knot  record,  entitled 
run  24.  The  response  points  contained  in  this  data  set  are  also  displayed  in  Table  7.1. 
Horizontal  tail  span  is  defined  from  the  right  tip.  Vertical  tail  span  begins  at  the  base  of 
the  tail. 


Solution  Parameters 


Response  Data 

Force  Application  Points 

Run  24:  Accelerometers 

Aircraft  Location 

%  Span 

NASTRAN 

Grid 

Nose  Cone  Lateral 

Horizontal  Tail 
C/4 

12.38 

48008 

Nose  Cone  Vertical 

28.39 

Pilot  Heel  Lateral  ■ 

45.84 

| 

Pilot  Floor  Lateral 

58.97 

44008  | 

Pilot  Floor  Vertical 

Pilot  Overhead  Lateral 

Copilot  Overhead  Lateral 

Vertical  Tail 

C/4 

80.66 

31844 

Vertical  Pylon  Top  Lateral 

31842 

Hor  Stab  Left  Tip  Long 

49.45 

31840 

Hor  Stab  Right  Tip  Long 

33.85 

31838 

Hor  Stab  Left  Tip  Vertical 

18.25 

31836 

Hor  Stab  Right  Tip  Vertical 

6.81 

30722 

Table  7. 1 :  Response  Data  and  Force  Application  Points. 


The  first  solution  set  used  two  forces,  one  on  the  right  horizontal  stabilizer  at  mid 
semi-span  and  the  other  at  mid-span  of  the  vertical  tail.  Then,  additional  unknown  forces 
were  included,  with  different  response  data.  For  several  frequency  bands,  the  condition 
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number  of  the  FRF  matrix,  Hxy,  became  exceedingly  high.  This  occurred  when  the  FRFs 
were  very  small  or  FRFs  from  different  points  were  of  equal  amplitude.  In  these  cases,  the 
precision  of  MATLAB®  could  not  distinguish  between  numbers  and  interpreted  the  two 
rows  of  the  matrix  as  duplicates  of  each  other.  Therefore,  the  inverse  was  not  reliable. 

Using  between  two  and  twelve  forces,  the  response  locations  were  varied  to 
produced  different  solutions.  Table  7.2  through  7.8  display  the  results  for  seven  different 
loading  conditions.  In  all  these  cases,  every  response  location  that  run  24  contained  was 
used.  The  tables  present  the  RMS  for  each  force  within  the  frequency  bands  previously 
identified.  In  every  case,  the  energy  within  these  three  bins  represented  over  80%  of  the 
total  energy,  throughout  the  5-25  Hz  range,  of  the  force  at  that  node.  The  final  column 
indicates  the  RMS  of  the  total  force. 

The  RMS  values  were  determined  for  each  frequency  band  and  then  divided  by 
the  largest  RMS  value  of  total  force.  They  represent  the  percentage  of  the  total  force  that 
is  located  at  that  point  on  the  aircraft  and  within  the  frequency  range.  Similarly,  the  RMS 
of  total  force  was  normalized  to  this  same  value.  Hence,  the  RMS  of  total  force  that 
equals  1.0  represents  the  location  on  the  tail  that  had  the  largest  force  throughout  the  5-25 
Hz  band.  Normalizing  the  forces  was  done  to  aid  in  comparison  of  loads  between  cases. 
Additionally,  the  conditioning  of  the  matrices  produced  highly  inaccurate  results  within 
certain  frequency  bands.  In  these  cases,  the  solution  had  large  amplitude  PSDs  in 
frequency  ranges  containing  little  response  and  small  FRFs.  The  loads  in  nine  and  twelve 
force  cases,  Table  7.7  and  7.8,  are  highlighted  to  show  the  frequency  content  and 
magnitudes  of  the  computed  loads. 
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Forcing  Function  -  Normalized  (Ib/lb) 

Frequency  Bin 

RMS  of 
Total  Force 

Node 

5-8  hz 

10-13  hz 

22-25  hz 

48008 

0.34 

0.38 

0.21 

0.44 

31840 

0.47 

1.00 

0.52 

1.00 

Table  7.2:  RMS  of  Two  Force  Solution  Set. 


Forcing  Function  -  Normalized  (lb 

>/lb) 

Frequency  Bin 

RMS  of 
Total  Force 

Node 

5-8  hz 

10-13  hz 

22-25  hz 

42008 

0.36 

0.42 

0.11 

0.60 

46008 

0.39 

0.39 

0.20 

0.71 

31840 

0.46 

0.71 

0.30 

1.00 

Table  7.3:  RMS  of  Three  Force  Solution  Set. 


Forcing  Function  -  Normalized  (Ib/lb) 

Frequency  Bin 

RMS  of 

Node 

5-8  hz 

10-13  hz 

22-25  hz 

Total  Force 

43008 

0.14 

0.06 

0.12 

0.22 

46008 

0.15 

0.19 

0.09 

0.28 

31844 

0.51 

0.49 

0.36 

0.88 

31836 

0.76 

0.38 

0.26 

1.00 

Table  7.4:  RM 

IS  of  Four  Force  Solution  Set. 
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Forcing  Function  -  Normalized  (Ib/lb) 

Frequency  Bins 

RMS  of 
Total  Force 

Node 

5-8  hz 

10-13  hz 

22-25  hz 

43008 

0.22  1 

0.15 

0.07 

0.30 

45108 

0.53 

0.28 

0.06 

0.63 

|  46008 

0.30 

0.16 

0.04 

0.36 

|  31844 

0.50 

0.16 

0.10 

0.60 

31836 

0.89 

0.21 

0.13 

1.00 

Table  7.5:  RMS  of  Five  Force  Solution  Set. 


Forcing  Function  -  Normalized  (Ib/lb) 

Frequency  Bins 

RMS  of 

Node 

5-8  hz 

10-1 3  hz 

22-25  hz 

Total  Force 

42008 

0.03 

0.01 

0.00 

0.03 

44008 

0.07 

0.02 

0.01 

0.07 

46008 

0.03 

0.01 

0.00 

0.03 

31844 

0.48 

0.11 

0.05 

0.47 

31840 

1.00 

0.27 

0.09 

1.00 

31836 

0.53 

0.17 

0.05 

0.54 

- 

fable  7.6:  RMS  of  Six  Force  Solution  Set. 

Forcing  Function  -  Normalized  (Ib/lb) 


Node 


Frequency  Bin 


5-8  hz 


10-13  hz 


22-25  hz 


RMS  of 
Total  Force 


46008 


0.00199 


0.00564 


0.00034 


0.00613 


45108 


0.01034 


0.00795 


0.00135 


0.01426 


44008 


0.01581 


0.01946 


0.00173 


0.02596 


43008 


0.00661 


0.01301 


0.00036 


0.01471 


31844 


0.19509 


0.15047 


0.00907 


0.25202 


31842 


0.40794 


0.32484 


0.01723 


0.53204 


31838 


0.58917! 


0.56503 


0.01648 


0.83418 


31836 


0.60257 


0.75681 


0.02316 


1.00000 


30722 


0.23072 


0.36292 


0.01422 


0.44914 


Table  7.7:  RMS  of  Nine  Force  Solution  Set. 
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Forcing  Function  -  Normalized  (Ib/lb) 

Frequency  bin 

RMS  of 
Total  Force 

Node 

5-8  hz 

10-13  hz 

22-25  hz 

48008 

0.00027 

0.00777 

0.00005 

0.00818 

46008 

0.00053 

0.02362 

0.00015  ! 

0.02487 

45108 

0.00047 

0.02538 

0.00018 

0.02738 

44008 

0.00148 

0.01474 

0.00021 

0.01830 

43008 

0.00122 

0.01450 

0.00018 

0.01670 

42008 

0.00028 

0.00363 

0.00004 

0.00415 

31844 

0.03813 

0.09873 

0.00017 

0.10854 

31842 

If  0.20775 

0.00037 

0.53079 

31840 

0.86229 

0.00037 

1.00000 

31838 

■fM 

■Si 

0.00030 

0.86572 

31836 

0.17668 

0.00009 

0.36249 

30722 

0.03467 

0.07152 

0.00010 

0.08449 

Table  7.8:  RMS  of  Twelve  Force  Solution  Set. 


Figure  7.1  through  7.4  show  the  PSDs  of  loads  for  the  twelve  force  solution  set. 
The  loads  on  the  horizontal  tail  are  within  frequency  bins  which  have  little  response.  The 
magnitudes  are  much  smaller  than  forces  on  the  vertical  and  associated  with  poorly 
conditioned  Hxy  matrices.  Under  the  PSD  plots  is  the  phase  angle  between  point  31840, 
where  the  force  of  largest  magnitude  is  located,  and  selected  points.  They  were 
determined  from  CSDs  which  capture  the  phasing  information.  Bendat  and  Piersol  [Ref. 
8:p.  125]  provide  a  method  of  determining  the  phase  angle  from  cross-spectral  density 
functions.  Although  the  plots  are  difficult  to  interpret,  if  viewed  at  just  the  modal 
frequencies  earlier  identified,  6.0  and  1 1.6  Hz,  the  phase  relationship  can  be  seen.  This  is 
available  on  the  vertical  tail  where  forces  at  these  frequencies  are  located.  The  horizontal 
tail  has  no  forces  at  these  frequencies. 
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zq/7qi  -  asd  pszhbiujon 


Phase  Angle  31 840  -  451 08 


-(80  1 ‘  '  i Li i i i i i i i 1 1 1 1 1 > * i * 1 

5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25 


Frequency  -  Hertz 

Figure  7. 1 :  PSD  of  Forces  on  the  Right  Horizontal  Tail. 
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Normalized  PSD  -  lb  /hz 


x  10 


-4 


PSD  of  Forces  -  Left  Horizontal  Stab 


5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25 

Frequency  -  Hertz 

Figure  7.2:  PSD  of  Forces  on  the  Left  Horizontal  Tail. 
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degrees  Normalized  PSD  -  lb  /hz 
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Normalized  PSD  -  lb  /hz 


PSD  of  Forces  -  Vertical  Stab  (0%-50%  Span) 
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Frequency  -  Hertz 


Figure  7.4:  PSD  of  Forces  on  the  Vertical  Tail  -  0%  to  50%  Span. 
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D.  ANALYSIS 


Beginning  with  the  initial  load  set  and  continuing  through  twelve  forces,  the 
computations  place  the  loads  of  largest  magnitude  on  the  vertical  tail.  Initially,  the  10-13 
Hz  frequency  bin  had  the  greatest  energy.  As  more  and  more  forces  were  added,  the 
magnitude  of  forces  in  the  5-8  Hz  range  grew  tremendously.  This  coincided  with  a  poorly 
conditioned  FRF  (Hxy)  matrix  at  5  Hz.  Even  when  eliminating  the  contribution  to  RMS 
from  this  frequency,  forces  from  the  base  of  the  vertical  tail  to  mid-span,  at  both  5-8  and 
10-13  Hz,  play  the  dominant  role  in  producing  the  response. 

In  going  from  nine  to  twelve  forces,  the  10-13  Hz  frequency  bin,  again,  became  the 
dominant  source  of  energy.  The  problems  with  the  conditioning  of  the  FRF  matrix, 
discussed  above,  were  again  present  in  the  final  solutions.  The  point  where  the  largest 
forces  were  placed  was  grid  31840.  Figures  7.5  and  7.6  display  the  distribution  of  forces 
along  the  horizontal  and  vertical  tail.  It  can  be  seen  that  more  energy  was  placed  on  the 
lower  two-thirds  of  the  vertical  tail  in  both  the  5-8  and  10-13  Hz  frequency  band.  This  is 
consistent  with  the  results  from  the  nine  load  case,  shown  in  Table  7.7. 

The  results  also  match  the  information  from  the  analysis  of  the  frequency  response 
functions  and  pressure  transducers.  The  solution  places  loads  on  the  fuselage  at  points 
that  have  high  magnitude  FRFs  at  1 1 .6  Hz.  Pressure  transducers  confirm  broadband 
excitation,  including  this  frequency,  and  of  largest  amplitude  at  the  base  of  the  vertical  tail. 

Although  the  solution  indicates  a  force  on  the  vertical  tail  in  the  5-8  Hz  frequency 
bin,  pressure  transducers  do  not  support  this.  The  solution  also  has  no  forces  in  the  22-25 
Hz  range.  The  lack  of  NASTRAN  modes  in  this  bin  has  previously  been  discussed. 
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RMS  Force  -  Normalized  (lb/lb]  ^  RMS  Force  -  Normalized  (Ib/lb] 


Vertical  Tail  (c/4)  -  RMS  of  Forces 


6.81%  18.25%  33.85%  49.45%  65.05%  80.66% 

%  Span  -  Bottom  of  Stab  to  Top 


Figure  7.6:  Distribution  of  Forces,  12  Force  Solution  Set  -  Vertical  Tail. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  goal  of  this  work  was  to  quantify  the  magnitude  and  location  of  aerodynamic 
loads  contributing  to  the  random  excitation  of  the  empennage.  Then,  follow-on  efforts 
would  have  the  proper  focus  and  direction  in  modifying  the  structure.  This  was 
accomplished  by  first  conducting  an  extensive  analysis  of  flight  test  response  data  and 
pressure  transducers.  Due  to  the  nature  of  the  accelerometer  time  histories,  spectral 
analysis  was  used.  With  this  accomplished,  the  mode  shapes  and  frequency  response 
functions  of  the  aircraft  were  determined.  A  NASTRAN  model,  modified  to  replicate  the 
test  flight  aircraft,  provided  this  information.  Finally,  a  MATL AB®  program  was  written 
to  compute  the  spectral  content  of  the  aerodynamic  loads  producing  the  random 
excitation. 

Each  portion  of  the  analysis  provided  insight  into  the  nature  of  the  response, 
locations  of  excitations,  and  the  method  in  which  the  structure  transfers  input  to  output. 
From  all  of  these  elements,  combined,  conclusions  were  drawn. 

The  dominant  response  is  lateral  and  occurs  in  three  frequency  bands  or  bins.  The 
first  bin  is  5-8  Hz,  then  10-13  Hz,  and  finally  22-25  Hz.  Each  accelerometer  has  its 
maximum  value  within  one  of  these  bins.  The  top  of  the  vertical  pylon  has  maximum 
response  at  6.00  Hz.  For  the  nose  cone,  pilot’s  station,  and  horizontal  tail,  10-13  Hz 
excitation  is  the  largest.  Finally,  the  22-25  Hz  response  is  restricted  primarily  to  the  CPG 
station  and  vertical  acceleration  of  the  horizontal  tail. 
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Within  the  first  two  bins,  there  are  mode  shapes  that  correspond  with  the  motion 
seen  in  the  output  of  the  22  accelerometers.  The  aft  fuselage  lateral  torsion  mode  at  6.00 
Hz  and  horizontal  stabilizer  yaw  mode  at  1 1.6  were  identified  not  only  through 
NASTRAN  but  also  during  actual  frequency  response  hover  testing.  No  mode  was 
identified  or  significant  frequency  response  function  found  within  the  22-25  Hz  bin.  The 
horizontal  stabilizer  vertical  bending  mode  at  26.02  Hz  matches  the  response  seen  in  the 
22-25  Hz  bin.  The  actual  natural  frequency  of  this  mode  was  not  identified  during  flight 
test;  however,  in  previous  NATRAN  models  it  is  at  22.06  Hz. 

From  the  forces  computed  through  the  MATLAB®  program,  good  correlation 
existed  with  the  previous  analysis.  Numerous  loading  combinations  were  applied  with 
maximum  loads  continuously  placed  on  the  vertical  stabilizer  from  the  base  to  50%  span. 
The  10-13  Hz  frequency  bin  contained  the  largest  force.  This  agreed  with  data  from 
differential  pressure  transducers  which  showed  the  highest  broadband  pressures  at  the  base 
of  the  vertical  tail. 

From  this,  it  can  be  concluded  that  forces  within  the  10-13  Hz  bandwidth 
impinging  on  the  lower  two-thirds  of  the  vertical  tail  play  the  most  significant  role  in 
creating  the  random  excitation.  This  is  supported  by  the  differential  pressure  transducers. 
Solutions  want  to  place  5-8  Hz  forces  in  this  same  area.  This  is  not  supported  by  pressure 
transducers  and  leads  to  the  conclusion  that  main  rotor  IP  is  the  primary  source  of  6  Hz 
excitation.  Because  the  22-25  Hz  excitation  is  so  localized  and  no  FRFs  of  magnitude  are 
in  this  region,  the  least  squares  solution  provides  no  contributing  loads.  If  the  horizontal 
stabilizer  vertical  bending  mode  is  not  at  26.02  Hz  but  lower,  the  broadband  excitation 
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seen  in  the  differential  pressure  transducers  can  excite  this  mode  and  produce  the  recorded 
response. 

B.  RECOMMENDATIONS 

1 .  Correlate  with  Wind  T unnel  Testing 

Boeing  and  Sikorsky  have  conducted  two  wind  tunnel  tests.  The  purpose  was  to 
quantify  the  flow  impinging  on  the  tail  and  identify  changes  to  the  surface  of  the  aircraft, 
such  as  strakes,  that  would  modify  the  flow  and  reduce  excitation  levels.  The  results  of 
this  analysis  should  be  compared  to  data  from  the  wind  tunnel.  This  could  confirm  the 
conclusions  of  this  study  about  the  nature  of  the  aerodynamic  loads  and  their  primary 
locations  on  the  fuselage. 

2.  Parameter  Identification  and  Correlation  with  NASTRAN 

During  hover  frequency  response  testing,  actual  natural  frequencies  were  identified 
for  three  modes.  As  additional  test  flight  hours  become  available,  more  parameter 
identification  can  occur.  This  will  better  identify  the  exact  frequency  of  the  horizontal 
stabilizer  vertical  bending  mode,  as  well  as  clarify  other  modes.  In  addition,  the  change  in 
natural  frequencies  due  to  adding  stiffening  brackets  and  tip  weights  to  the  horizontal 
stabilizer  can  be  assessed.  With  this  information,  a  detailed  analysis  of  the  NASTRAN 
model  can  be  conducted  with  modification  of  natural  frequencies  and  mode  shapes  to 
replicate  the  actual  aircraft. 
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3.  Additional  Airspeed  Analysis 

The  computed  solution  using  the  method  of  Bendat  and  Piersol  was  done  for  only 
100  knot  records.  After  receiving  additional  airspeed  records  from  Sikorsky,  the  PSDs 
and  CSDs  were  computed;  however,  no  complete  solutions  were  obtained  for  the  spectral 
content  of  forces.  Calculating  solutions  throughout  the  range  of  airspeeds  can  provide 
further  confirm  and  clarity  as  to  when  the  random  excitation  begins  and  if  the  6  Hz 
response  is  due  to  main  rotor  IP. 

4.  Structural  Optimization 

Ultimately,  the  goal  is  to  eliminate  or  reduce  the  response  to  acceptable  levels. 
Using  the  conclusions  from  this  analysis  and  confirmation  from  the  sources  discussed 
above,  the  structure  can  be  altered  to  reduce  the  frequency  response  functions  within  the 
excitation  bandwidth.  This  can  be  accomplished  using  several  different  methods,  or, 
optimally,  a  combination  of  all. 

First,  using  the  NASTRAN  model,  structural  modifications  can  be  done  that  move 
the  natural  frequency,  1 1 .6  Hz,  of  the  horizontal  stabilizer  yaw  mode.  The  focus  of  the 
changes  should  occur  on  the  lower  two-thirds  of  the  vertical  tail.  Extensive  work  was 
done  by  Tobin  [Ref.  12]  and  Shoop  [Ref.  13]  to  increase  torsional  stiffness  of  the  aft 
tailcone.  The  work  attempted  to  stiffen  the  tail  and  move  the  aft  fuselage  lateral  torsional 
mode  away  from  MR  IP.  Their  analysis  and  strain  energy  diagrams  provide  insight  to 
methods  of  altering  the  structure,  specifically  for  the  1 1.6  Hz  mode. 

Second,  changes  are  proposed  to  the  design  of  the  tail  fold,  spars  of  the  vertical 
tail,  and  the  tail  fan  gear  box  struts.  Additionally,  a  trade  study  is  being  conducted  by  the 
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Comanche  Program  Management  Office  (PMO  -  Comanche)  on  the  mounting  of  the  tail 
landing  gear  bay.  As  these  designs  are  finalized,  the  dynamics  of  the  tail  will  change.  The 
quantification  of  the  flow  determined  through  this  analysis  should  aid  in  creating  a 
dynamically  “tuned”  design. 

Finally,  modifying  the  flow  around  the  main  rotor  hub  and  pylon  has  been 
investigated  through  the  wind  tunnel  experiments  previously  discussed.  Modifications  to 
the  structure  must  be  done  in  conjunction  with  changes  to  the  flow.  Only  through  this 
integrated  approach  can  an  acceptable  solution  be  obtained.  In  this  manner,  both  the  flow 
and  structural  response  can  be  optimized  for  maximum  reduction  of  the  airframe  vibratoiy 
response. 
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APPENDIX  A:  TEST  FLIGHT  PARAMETERS 


Parameter  y- 

Affected  Components 

.Units 

HBN0RML1 

Horz  Stabilizer  Flatwise 

Bending,  LHS  (near  Root) 

Horizontal  Stabilizer 

in-lb 

HBNORMR1 

Horz  Stabilizer  Flatwise 

Bending,  RHS  (near  Root) 

Horizontal  Stabilizer 

in-lb 

HBNORML2 

Horz  Stabilizer  Flatwise 

Bending,  LHS  (outboard) 

Horizontal  Stabilizer 

N/A 

HBNORMR2 

Horz  Stabilizer  Flatwise 

Bending,  RHS  (outboard) 

Horizontal  Stabilizer 

N/A 

HSTBFWDB 

Horz  Stab  Fwd  Folding  Joint 
Bending  LHS 

Horizontal  Stabilizer 

in-lb 

HSVAFTL 

Horz  Stab  Aft  LHS  Spar  Shear 
Force  -  lb 

Horizontal  Stabilizer 

N/A 

HSVAFTR 

Horz  Stab  Aft  RHS  Spar  Shear 
Force  -  lb 

Horizontal  Stabilizer 

N/A 

HSVFWDL 

Horz  Stab  Forward  LHS  Spar 
Shear  Force  -  lb 

Horizontal  Stabilizer 

N/A 

HSVFWDR 

Horz  Stab  Forward  RHS  Spar 
Shear  Force  -  lb 

Horizontal  Stabilizer 

N/A 

VSFFTS1 

Vert  Stab  Attach  Fwd  Ftg  Outer 
Cap  LHS 

Vert  /  Shroud 

u‘7” 

VSFFTS2 

Vert  Stab  Attach  Fwd  Ftg  Outer 
Cap  RHS 

Vert  /  Shroud 

u‘7” 

VSAFTS1 

Vert  Stab  Attach  Aft  Ftg  Outer 
Cap  LHS 

Vert  /  Shroud 

u‘7” 

VSAFTS2 

Vert  Stab  Attach  Aft  Ftg  Outer 
Cap  RHS 

Vert  /  Shroud 

u ar 

Table  A.  1 :  Test  Flight  Parameters. 
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|  Mnemonic 

Parameter 

Affected  Components 

Units 

BVSRTMX 

Vertical  Stab  Root  Roll 
Moment 

Vert  /  Shroud 

DHSDRAG 

Horizontal  Stab  Drag  -  DP 

Horizontal  Stab 

EB 

DHSTABMX 

Net  Stab  Roll  Moment  -DP 

Horizontal  Stab 

DHYAWMZ 

Horizontal  to  Vertical  Stab 

Yaw  Moment  -  DP 

Horizontal  Stab 

DHSTABVZ 

Net  Stabilizer  Lift  (lb)  -  DP 

Horizontal  Stab 

N/A 

DVSAFTLH 

Derived  Vert  Load  Aft  LHS 
Lug  -  DP 

Vert  /  Shroud 

DVSAFTRH 

Derived  Vert  Load  Aft  RHS 
Lug  -  DP 

Vert  /  Shroud 

lb 

DVSFWDLH 

Derived  Vert  Load  FWD  LHS 
Lug  -  DP 

Vert  /  Shroud 

lb 

DVSFWDRH 

Derived  Vert  Load  FWD  LHS 
Lug  -  DP 

Vert  /  Shroud 

lb 

DVSRTMX 

Vertical  Stab  Root  Roll  MMT 
-DP 

Vertical  Stab 
Attachment 

N/A 

DVSRTMY 

Vertical  Stab  Root  Pitch 

MMT -DP 

Vertical  Stab 

Attachment 

N/A 

DVSTABVY 

Derived  Vertical  Stab  Lateral 
Force  -  DP 

Vert  /  Shroud 

N/A 

HBEDGEL1 

Horz  Stab  Edgewise  Bending, 
LHS  (near  root) 

Horizontal  Stab 

N/A 

HBEDGER1 

Horz  Stab  Edgewise  Bending, 
RHS  (near  root) 

Horizontal  Stab 

N/A 

Table  A.  2:  Test  Flight  Parameters. 
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APPENDIX  B:  PSD  AND  CSD  PROGRAMS 


%  psd_flt__data.m 

%  Program  reads  the  input  flight  data  for  a  run  and  the  computes  the 
%  PSD  for  each  accelerometer  time  history 

%  The  input  record  is  an  18,000  X  n  matrix  with  each  column 
%  representing  another  transducer 


%  Input  the  data  set  you  want  to  read,  sampling  rate,  data  set  you 
%  want  plotted 

clear 

run_labels__file 

dt  =  1/300; 
string  = [ 1 runl40 T ] ; 
speed  =[*140  Knots1]; 
eval ( [ *  load  ' , string, ';']); 
i  =  0; 

for  col  =  4:32 
i  =  i  +  1; 
n  =  1024; 

noverlap  =  round (. 67*n) ; 

eval([*x  =  1 , string, *_data (:,  col)  ;’])  ; 

%  Using  the  SPEC  command 

%  [Pxx,F]  =  spec (x, n, noverlap, hanning (n) , 1/dt) ; 

%  Using  the  PSD  command 

[Pxx,F]  =  psd (x, n, 1/dt, hanning (n) , noverlap) ; 

%  Normalize  by  multiplying  by  2/frequency  sampling  rate 
%  To  get  the  rms  multiply  by  sampling  rate/nfft  then  add 
%  This  is  from  Pareseval's  theorem 

Pxx  =  (2/ (1/dt) )  *  Pxx; 

%  Another  way  to  normalize  so  that  the  rms  is  just  amplitudes  add 
%  together 

%  With  this  method,  the  value  from  different  length  fftTs  is  the  same 
%  Pxx  =  (2/ (1/dt))  *  ( 300/n)  *  Pxx; 
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%  Plotting  Function 


plot_max  =  round(25/ ( (1/dt) /n) ) ; 

ymax  =  max (Pxx ( 10 : plot_max, 1) ) ; 

figure (i) 

plot (F, Pxx ( : , 1) ) ; 

grid; 

axis ( [0  25  0  ymax] ) ; 

set (gca, ’xtick' , [0 : 1:25] )  ; 

eval ( [ ' title ( ' ' PSD  deblank (run_labels (col, num2str(n),'  Point 
FFT  ”)']); 

xlabel (' Frequency  -  Hz'); 
ylabel ( 'g'sA2/Hz' ) ; 

eval ( [ 'text (5, . 9*ymax, ' ' Run  ' , speed, . . . 

',  Hanning  Window,  67%  Overlap'')']); 

%  Save  the  PSD  data  in  another  matrix  for  use  in  solving  for  forcing 
functions 

%  Name  of  the  matrix  saved  is  psd_run24  (mat  format)  and  has  the 
frequency  and 

%  the  psd  for  all  accelerometer  data  we  had  for  run24. 


eval ([ ’psd_' , string, '(:,  i)  =  Pxx ( : ,  1) ; '  ]  )  ; 
end 

%save  frequency  as  a  separate  variable 
psd_labels  =  run_labels (4 : size (run_labels, 1) , : ) ; 

eval ([ 'save  psd_', string, '  psd_' , string, '  F  psd_labels ' ] ) ; 
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%  csd_flt_data .m 

%  Program  reads  the  input  flight  data  for  a  run  and  the  computes  the 
%  CSD  (cross  spectral  density)  between  output  accelerometer  time 
%  history 

%  The  input  record  is  an  18,000  X  n  matrix  with  each  column 
%  representing  another  transducer 


%  Input  the  data  set  you  want  to  read,  sampling  rate,  data  set  you 
%  want  plotted 

dt  =  1/300; 
load  run24 
k  =  0; 

for  col  =  5:5 
for  i  =  1:12 
k  =  k  +  1; 
ii  =  i  +  3; 
n  =  1024; 

noverlap  =  round (. 67*n) ; 
x  =  run24_data ( : ,  col) ; 
y  =  run24__data  ( : ,  ii)  ; 


%  Using  the  CSD  command 

[Pxy, F]  =  csd(x,y,n, 1/dt, hanning (n) , noverlap) ; 

%  Normalize  by  multiplying  by  2/frequency  sampling  rate 
%  To  get  the  rms  multiply  by  sampling  rate/nfft  then  add 
%  This  is  from  Pareseval's  theorem 

Pxy  =  (2/ (1/dt))  *  Pxy; 

%  Another  way  to  normalize  so  that  the  rms  is  just  amplitudes  added 
%  together 

%  With  this  method,  the  value  from  different  length  fft's  is  the  same 

%Pxy  =  (2/ (1/dt))*  ( 300/n)  *  Pxy; 

%  Plotting  Function 

plotjnin  =  round ( 5/ (( 1/dt ) /n) ) ; 

plot_max  =  round (25/ ( (1/dt) /n) ) ; 

ymax  =  max ( abs ( Pxy (plot_min : plotjmax, 1 ) ) ) ; 

figure (k) 

subplot (2,1,1)  ; 

plot ( F, abs ( Pxy ( : , 1 ) ) ) ; 

grid; 

axis  ([5  25  0  ymax]); 
set (gca, 'xtick' , [5:1:25] ) ; 

eval ( [ 'title ( 1 ' CSD  ' , deblank ( run24_labels (col,  :  ) ) ,  1  -  ' , 
deblank (run24_labels  (ii, : ) ) ,  1  1  1 )  1  ]  )  ; 

%xlabel ( ' Frequency  -  Hz'); 
ylabel ( 'g'sA2/Hz  (Amplitude) '); 
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subplot  (2, 1,2) 

plot ( F, 180/pi  *  angle (Pxy (:, 1) )) ; 

angle_max  =  max  (180/pi  *  angle  (  Pxy  (plot_min:plot__max,  1 )))  ; 

axis  (  [5  25  -180  180] )  ; 

set (gca, *  xtick 1  ,  [5:1:25]); 

set (gca, 1 ytick1 ,[-180:60:180]); 

grid; 

xlabel ( 1  Frequency  -  Hz’); 
ylabel ( 1  Phase  Angle  -  degrees1) 

eval ( [ 1  title ( 1 1  Run  1 , num2str (run24_data (1,1)),... 

f,  Hanning  Window,  67%  Overlap,  1 , num2str (n)  ,  1  Point  FFT 1 1 )  1 ] )  ; 
%eval ( [ 1  text (5,  . 8 9*angle_max, 1 1  Phase  Angle 
1 ,  num2str ( run24_data (1,1)),... 

%  1 ,  Hanning  Window,  67%  Overlap,  1 , num2str (n)  ,  1  Point  FFT11)1]) 

%  Save  the  CSD  data  in  another  matrix  for  use  in  solving  for  forcing 
%  functions 

%  Name  of  the  matrix  saved  is  psd_run24  (mat  format)  and  has  the 
%  frequency  and  the  psd  for  all  accelerometer  data  we  had  for  run24. 

eval  (  [  1  csd_run24__* ,  num2str  (col-3 )  ,  1  ( : ,  i )  =  Pxy  (:,1);*]); 

end 


%eval ( [ 1  save  psd_run24  csd_run24_’ , num2str (col-3 ),  1  F  -append;1]); 
end 
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APPENDIX  C:  MODE  SHAPES 


Mode  shapes  for  selected  modes  are  presented  on  the  next  eleven  pages. 


123 


Aft  Fuselage  Lateral  Torsional  Mode,  6.00  Hertz 
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.004115 


Horizontal  Stab  Roll  Mode,  9.50  Hertz 


.003876  I  _  II  . . I  . I  _ .003876 


Horizontal  Stab  Yaw  Mode,  1 1 .60  Hertz 


.002863 


2nd  Fuselage  Vertical  Bending  Mode,  14.52  Hertz 


2nd  Fuselage  Lateral  Mode,  16.07  Hertz 


.006288 


Fuselage  Mode,  24.04  Hertz 


.001806 


Horizontal  Stab  Vertical  Bending  Mode,  26.02  Hertz 


APPENDIX  D:  ACCELEROMETER  LOCATIONS  &  NASTRAN  GRIDS 


RAH66- MODEL  J  yrahj194.dat  4/17/98  sp 

COMPARISON  OF  KEY  FLIGHT  TEST  LOCATIONS  WITH  SUGGESTED  NASTRAN  APPROXIMATE  LOCATIONS 


COPILOT  HEEL 


5  COPILOT  FLOOR 


6  PILOT  OVERHEAD 


COPILOT  OVERHEAD 


8  VERTICAL  PYLON  TOP 


HOR.  STABILIZER  LEFT  TIP  VSTABTLT:D:300  17670.00  -1100.00  5463.00 
LOSTABTL:D:3CO 


11  HOR.  STABILIZER  VERTICAL  HBNORML1:D:300l 
BEND.  LT  SIDE 


FLIGHT  TEST 
Approximate  Location 

Mnemonic  STA  BL  WL  ST, 

mm  mm _ mm 


VNOSE:D:300  5286.00  0.00  2720.00 

LNOSE:D:300 
LONOSE:D:300 

LPHEELD300  6452.00  150.00  2720.00 

VPHEELD300 


LPF:D300  7450.00  0.00  2720.00 

LOPF:D300 
VPRD300 


VCOPHEEL:D:300 1  7740.00  150.00  2980.00 


VCOPF:D:300  8900.00  0.00  3130.00 


LPOVRHD:D:300  7420.00  0.00  4100.00 


LCPOVRH D: D: 300 1  8795.00  200.00  4200.00 


LVPYLNT:D300  17737.00  185.00  5360.00 


10  HOR.  STABILIZER  RIGHT  TIP  VSTABTRT:D:300  17670.00  1580.00  5463.00 

LOSTABTR:  D:300 
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FLIGHT  TEST 

NASTRAN 

NASTRAN 

Approximate  Location 

Location 

GRID 

TEM 

LOCATION 

Mnemonic 

STA  BL  WL 

STA  BL  WL 

POINT 

No. 

mm  mm  mm 

mm  mm  mm 

12 

HOR.  STABILIZER  VERTICAL 

HBNORML2:D:300 

-660  00 

17670.10  -64836  5463.39 

42108 

BEND.  LT  SIDE 

17666  70  -648  36  5374  97 

42101 

17982  30  -648.36  5444.27 

m 

17979  90  -648  36  5367.95 

13 

HOR.  STABILIZER  VERTICAL 

HBNORMR1  :D:300 

380.00 

17670.10  341.60  5463.39 

45308 

BEND.  RT  SIDE 

17666.70  341.60  5375  00 

45301 

17982  40  415.31  5444.25 

45305 

17980.00  415.31  5367.95 

45304 

14 

HOR.  STABILIZER  VERTICAL 

HBNORMR2:D:300 

770  00 

17670.10  656.83  5463  39 

46008 

BEND.  LT  SIDE 

17666  70  656.83  5374.97 

46001 

17982.60  656.83  5444.24 

46005 

17979  90  656.83  5367.95 

46004 

15 

HOR.  STABILIZER  EDGE 

HBEDGEL1:D:300 

-240  00 

17670.10  -281.41  5463  39 

43108 

BEND.  LT  SIDE 

17666.70  -281.41  5374.97 

43101 

17982.30  -281.41  5444.27 

43105 

17979  90  -281.41  5367.95 

43104 

16 

HOR.  STABILIZER  EDGE 

HBEDGER1  :D:300 

350  00 

17670.10  341.60  5463.39 

45308 

BEND.  RT  SIDE 

17666.70  341.60  5375.00 

45301 

17982  40  415.31  5444.25 

45305 

17980.00  415.31  5367.95 

45304 

17 

LT  AFT  SHROUD/V.TAIL 

VSAFTS1:D:300 

17002.00  -294  50  4400.00 

21326 

INTERFACE  FTTING 

18 

RT  AFT  SHROUD/V.TAIL 

VSAFTS2:D:300 

17002  00  -176.00  4400.00 

21327 

INTERFACE  FTTING 

19 

LT  FWD  SHROUD/V.TAIL 

VSFFTS1:D:300 

16635  00  -294  50  4400  00 

24621 

INTERFACE  FTTING 

20 

LT  AFT  SHROUD/V.TAIL 

VSFFTS2:D:300 

16635  00  -176  00  4400.00 

24620 

INTERFACE  FTTING 

21 

FWD  H.  STAB/V.TA1L  ■ 

DH  STAB  MX:  D300 

17730  00  62.85  5404. 00 ' 

32391 

FITTING  CL 

DHSTABMZ:D300 

17730.00  151.15  5404.001 

32394 

DHSTABVZ:  D300 

17737.10  104  55  5360  00 

32336 

Table  D.2:  Accelerometer  -  NASTRAN  Locations. 
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APPENDIX  E:  LOAD  CALCULATIONS  PROGRAM 


%  f  f_solve_nonsym_inter ,m 

%  Program  that  solves  for  the  forcing  functions,  in  PSD  by  using  the 

%  techniques  discussed  by  Bendat  and  Pierson  for  multiple  input 

%  single  output 

%  Gxx  =  (Hxy* 1 ) ^-1  *  Gyy  *  (Hxy) A-1 


%  First  must  identify  the  location  (nodes)  of  the  forcing  function 

%  the  location  of  the  response  (nodes),  and  the  direction  of  the 

%  response 

%  The  directions  are  somewhat  backwards  -  due  to  the  way  I  had  Nastran 
%  write  out  the  punch  files  (  which  became  .mat  files) 


%  x  =  3  (longitudinal) 

%  y  =  2  (lateral) 

%  z  =  1  (vertical) 


%  The  response  data  uses  the  following  script: 

%  First  is  node  in  Nastran  Mode.l,  next  is  direction  and  last  is  the 
%  pointer  into  the  psd  matrix  (which  column) 


ncl  = 

[61641 

2 

l]  ; 

phi  = 

[6505 

2 

2] ; 

pfl  = 

[7246 

2 

3] ; 

pol  = 

[7514 

2 

4]  ; 

cpol  = 

[8614 

2 

5]  ; 

ncv  = 

[61641 

1 

6]  ; 

pfv  = 

[7246 

1 

Tit 

vptl  = 

[32341 

2 

8]; 

hrtv  = 

[48108 

1 

•  9]; 

hltv  = 

[41108 

1 

10]  ; 

hrtlo  = 

[48108 

3 

11]  ; 

hltlo  = 

[41108 

3 

12]  ; 

%  Nose  Cone  Lateral 
%  Pilot  Heel  Lateral 
%  Pilot  Floor  Lateral 
%  Pilot  Overhead  Lat 
%  CPG  Overhead  Lat 
%  Nose  Cone  Vertical 
%  Pit  Floor  Vertical 
%  Vert  Pylon  Top  Lat 
%  HorStab  Rt  Tip  Vert 
%  HorStab  Lt  Tip  Vert 

%  HorStab  Rt  Tip  Long 

%  HorStab  Lt  Tip  Long 


%  Use  the  following  format  to  enter  the  input-ouput  locations  and 
%  specify  which  run  set 

run  =  24; 

%no_fft  =  4096; 

%  force  =  [32311  46008  ] 

%  response  =  [  pol 
%  phi] ; 

force  =  [  42008  44008  46008  31844  31840  31836  ]; 
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dP  dP  dP  dP  dP 


response  =  [  phi 
pfv 
cpol 
ncl 
ncv 
vptl 
hrtv 
hltv 
hrtlo 
hltlo  ] ; 

%  Assumptions  is  that  loads  are  correlation  -  are  the  loads  correlated? 

%  If  the  loads  are  completely  uncorrelated  then  set  corr  =  0 

%  If  it  is  something  inbetween  then  define  the  matrix  relate 

force_corr  =  1; 
response_corr  =  1; 

%  Use  the  following  format  for  relate  -  forcing  function (load)  to  load 
(which  are  correlated) 
force_relate  =[110 
110 

001]  This  tells  us  force  1  is  related  to  force  2, 
all  others  are  not 

%  force_relate  =[10 

%  0  1]; 

%  Load  the  frequency  response  function  matrix 

%  Counter  to  use  for  matrix  index 

k  =  0; 


[no_row, force_num]  =  size(force); 

[ response^num, no__col ]  =  size  ( response)  ; 


for  f  =  l:force__num 

%  frf  file  must  have  the  name  format  frfnode 

for  r  =  1 : response_num 
k  =  k+1; 

eval(['load  f rf *, num2str ( force ( 1, f) ) , ' ;  '  ]) 

eval ( [ 1 f rf_sel  =  f rf_’ , num2str ( response ( r, 1 )),  1  ;  *]) 

y  =  [ ( f rf_sel ( : , 1 )  +  f rf_sel ( : , 2 )  .*i)  ( frf_sel ( : , 3)  + 

frf_sel(:,4)  .*i)  (frf_sel ( : ,5)  +  frfjsel ( : , 6)  .*i)]; 

H ( : , k)  =  y ( : f response ( r , 2 ) ) ; 
end 


end 
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%  H  (frf  matrix  has  now  been  loaded  in  the  following  format;  from  our 
%  example  nodes  - 

%  H  =  [Hll  H12  H13  H21  H22  H23  H31  H32  H33] ; 

% 

%  Because  of  the  weird  way  I  saved  the  pch  files,  the  vector  {time}  is 
%  actually  frequency 

%  Now  we  need  to  load  the  response  data  from  the  PSD  and  CSD  we  have 
%  done  before 


eval (['load  psd_run* ,num2str (run) , 1 ;  1 ]) 
if  exist  ('no_fft')  —  1 

eval ( [ ' load  psd^run' , num2str (run) ,  , num2str (no_f ft ) ,  1 ;  1 ] ) 

end 


k  =  0; 

for  f  =  1 : response_num 

for  r  =  1 : response_num 
k  =  k+1; 
if  f  ==  r 

G(:,k)  =  psd_run24 (:, response (r, 3) ) ; 
else 

eval  ( [ *G( : ,k)  =  ... 

csd_run24_'  ,num2str  (response  (f ,  3)  )  ,  '  ( : ,  response  (r,  3)  )  ;  '  ]  ) 

end 

end 

end 


%  Form  the  Hxy  and  Gyy  matrix  for  each  individual  frequency, 
%  solve  for  Gxx 
%  and  save  in  another  matrix 

kk  =  0; 

start_freq  =  5.50; 
end_freq  =  24.00; 

[f__dif f ,  f_index_start]  =  min  (abs  (start_f req  -  F)  )  ; 

[f__diff ,  f_index_end]  =  min  (abs  (end_f req  -  F)  )  ; 

for  ff  =  f_index__start  :  f_index_end 

G_count  =  0; 
k  =  0; 

kk  =  kk  +  1; 

for  1  =  1 : response_num 
for  m  =  1 : response_num 
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dP  dP  dP  dP  dP 


%  Load  the  Gyy  square  matrix 

Gyy(l,m)  =  G ( f f , G_count+m)  ; 
end 

G_count  =  G_count  +  response^num; 
end 

for  1  =  l:force_num 
for  m  =  1 : response_num 

Load  the  Hxy  square  matrix 

[min^dif , val_index]  =  min (abs (F ( f f )  -  time’)); 
Cubic  Spline  Interpolation 

Hxy(l,m)  =  interpl (time, H (:, k+m) , F ( ff) ,  1  spline  1 ) ; 
Linear  Interpolation 

Hxy (l,m)  =  interpl (time, H ( : , k+m)  ,  F ( f f ) , 1  spline 1 ) ; 
end 


k  =  k  +  response_num; 
end 

if  exist ( 1 f orce_relate 1 )  ==  1 
Hxy  =  force_relate  .  *  Hxy; 
end 

if  exist (' response_relate * )  ==  1 

Gyy  =  excite_relate  •*  Gyy; 
end 

if  response_corr  ==  0 

Gyy  =  eye (size (Gyy) )  .*  Gyy; 

end 


%  Now  we  have  to  solve  for  the  Gxx  matrix  via  the  technique 
%  discussed  above 

if  force^corr  ==  0  &  response_corr  ==  0 
Gyy  =  diag (Gyy)  ; 

Gxx  =  ( (abs (HxyT ) ) .A2)  \  Gyy; 
else 

Gxx  =  Hxy ’ \  Gyy  /Hxy; 
end 

%  Saving  the  output  data  Gxx  to  a  matrix  of  the  following  form 
%  Gxx  =  [Gil  G12  G13  G21  G22  G23  G31  G32  G33] 

%  Also  outputs  a  vector,  Freq,  that  has  the  frequencies 
%  each  row  represents 

[Gxx  row, Gxx_col]  =  size (Gxx); 


if  Gxx_col  ==  1 
G_output (kk, : )  =  Gxx1; 

else 
k  =  0; 
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for  11  =  l:Gxx_row 
for  mm  =  1 : Gxx__col 

G_output  (kk,  k-f-mm)  =  Gxx(ll,mm); 
end 

k  =  k+  force_num; 

end 

end 

frequency (kk)  =  F(ff); 
cond_Hxy ( kk )  =  cond ( Hxy ) ; 
rank_Hxy(kk)  =  rank (Hxy); 
end 

bigJKahuna  =  [frequency1  G_output  cond_HxyT  rank_Hxy’]; 
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